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ABSTRACT 

Context. This is the fourth paper in a series showing the results of planet population synthesis calculations. In Paper I, we presented 
our methods. In Paper II, we compared statistically the synthetic and the observed planetary population. Paper III addressed the 
influences of the stellar mass on the population. 

Aims. Our goal in this fourth paper is to systematically study the effects of important disk properties, namely disk metallicity, mass 
and lifetime on fundamental properties of planets like mass and semimajor axis. 

Methods. For a large number of protoplanetary disks which have properties following distributions derived from observations, we 
calculate a population of planets with our formation model. The model is based on the classical core accretion paradigm but includes 
self-consistently planet migration and disk evolution. 

Results. We find a very large number of correlations: Regarding the planetary initial mass function, metallicity, M disk and r disk have 
different roles: For high metallicities, giant planets are more frequent. For high Mdi sk , giant planets are more massive. For long t^, 
giant planets are both more frequent and massive. At low metallicities, very massive giant planets cannot form, but otherwise giant 
planet mass and metallicity are nearly uncorrelated. In contrast, (maximal) planet masses and disk gas masses are correlated. The 
formation of giant planets is possible for initial planetesimal surface densities S s of at least 6 g/cm 2 at 5.2 AU. The sweet spot for 
giant planet formation is at ~ 5 AU. In- and outside this distance, higher Ss are necessary. Low metallicities can be compensated 
by high M disk , and vice versa, but not ad infinitum. At low metallicities, giant planets only form outside the ice line, while at high 
metallicities, giant planet formation occurs throughout the disk. The extent of migration increases with M^sk and and usually 
decreases with metallicity. No clear correlation of metallicity and the semimajor axis distribution of giant planets exists because in 
low metallicity disks, planets start further out, but migrate more, while the contrary applies for high metallicities. The final semimajor 
axis distribution contains an imprint of the ice line. Close-in low mass planets have a lower mean metallicity than Hot Jupiters. The 
frequency of giant planets varies approximately as Mj L 2 k and TjL k . 

Conclusions. The properties of protoplanetary disks - the initial and boundary conditions for planet formation - are decisive for the 
properties of planets, and leave many imprints in the population. 

Key words. Stars: planetary systems - Stars: planetary systems: formation - Stars: planetary systems: protoplanetary disks - Planets 
and satellites: formation - Solar system: formation 



1. Introduction 

The number of known extrasolar planets has grown large enough 
to look at the statistical properties of the population as a whole, 
rather than at the properties of single objects, and to compare 
the actual population with a synthetic population obtained from 
a theoretical planet formation model. In this way all discov- 
ered planets serve (provided the detection bias is known) to con- 
strain the model, and to improve our understanding of planet for- 
mation. We used our extended core accretion formation model 
(Alibert et al. 2005a) to generate in a Monte Carlo way popu- 
lations of synthetic extrasolar planets (Mordasini et al. [2009a 
hereafter Paper I). Then, we compared the detectable synthetic 
planets with an observational comparison sample of actual exo- 
planets using statistical methods (Mordasini et al. |200 9b here- 
after Paper II). We found that we could reproduce in a statisti- 
cally significant way some of the most important properties of 
the observed extrasolar giant planets. Finally, in Alibert et al. 
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(I2010I I. hereafter Paper III, we discussed the influence of the stel- 
lar mass on the synthetic planetary population, studying for ex- 
ample the effect on the planetary initial mass function, the semi- 
major axis distribution or the "metallicity effect" (the increase 
of the detection probability of giant planets with metallicity). In 
this Paper IV, we focus back onto solar-like stars. 

Over the last years, observational considerable progress was 
achieved in the characterization of the end products of plane- 
tary formation process, i.e. the planet themselves. Progress has 
also been substantial in the characterization of the initial and 
boundary conditions for this process, i.e. the properties of pro- 
toplanetary disks. This was made possible to a large measure 
thanks to new observational facilities like Spitzer (e.g. Fang et al. 
2009 ). Observations of disks around young stars have provided 
us with knowledge of the distributions of disk masses (Beckwith 
& Sargent 1996, Andrews et al. 12009). disk sizes (McCaughrean 
& Odell 119961 Andr ews et al. 120101 ) and lifetimes (Haisch et 
al. 120011 Fedele et al. 120101) . Additionally, correlations between 
disk properties and stellar mass were discovered (Kennedy & 
Kenyon 2009i Mamajek 2009), which have important implica- 
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tions for the formation of planets around stars of different masses 
(Currie 2009 Paper III). As was shown in Paper I, the diversity 
of extrasolar planets is a direct consequence of the diversity of 
protoplanetary disks, which means that the properties of disks 
are critical in defining the outcome of planet formation taking 
place in such disks. 

In this article we focus on correlations between disk and 
planetary properties. We study systematically the influences of 
disk metallicity [Fe/H], disk (gas) mass and lifetime Tdisk 
on important, observable properties of synthetic planets by com- 
puting a large number of models of planet formation. 

1.1. Observed correlations 

From an observational point of view, a number of correlations 
were inferred in the past years, with various degrees of signif- 
icance (see Udry & Santos l2007l or Mayor et al. d201 11 1 for an 
overview). We address in this paper the following correlations: 

(1) The clearest correlation is the link between the stellar 
metallicity (to first order an indicator of the initial dust-to-gas- 
ratio in the disk, see Santos et al. 120031 ) and the likelihood of 
detecting a giant planet. This "metallicity effect" is observation- 
ally known for a long time and very well established for solar- 
like stars by numerous studies (e.g. Gonzalez 119971 Fischer & 
Valenti l20031 Udry & Santos l2007t . The "metallicity effect" was 
studied as an observational constraint in Paper II and is here fur- 
ther addressed in Sect.|42j|43]and|8]T] 

(2) Interestingly, lower mass Neptunian and Super-Earth 
planets seem in contrast not to be found preferentially around 
high [Fe/H] stars (Mayor & Udry [20081 Sousa et al. 120081 
Ghezzi et al. 1201 Oi l. We study this in Sect. [41] and [7] 

(3) Already much less compelling is a possible absence of 
very massive planets (in a mass range of about 5 to 20 Jupiter 
masses, M^) orbiting within a few AUs low metallicity solar- 
like stars as found in radial velocity (RV) searches. Such a possi- 
ble absence was noted by Udry et al. ([2002), Santos et al. (2003 ) 
and Fischer & Valenti (2005). This is discussed in Sect. |4.3.3] 

(4) In contrast, no secure correlations were found between 
stellar metallicity and planetary semimajor axis, at least among 
RV detections (Udry & San tos 120071 Valenti & Fischer 120081 
Ammler-von Eiff et al. 2009), even though possible correlations 
were discussed in the literature. Sozzetti (2004| for example 
studied whether stars with Pegasi planets have particularly high 
[Fe/H], even among planet hosts. Correlations between migra- 
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tion, semimajor axis and metallicity are addressed in Sect 
and [63] 

(5) Giant planets are more frequently found orbiting A stars 
than orbiting solar-like stars, and are also more massive (Lovis 
& Mayor 2007; Bowler et al. 120101) on average. In Paper III it 
was shown that this correlation can be best reproduced if one 
assumes a roughly linear scaling between disk (gas) mass and 
stellar mass. This implies, at least partially, a correlation between 
disk mass and planet mass (Sect. 4. 1.2 and 4.4 1. 

In addition, we address below several other correlations that 
appear in the models but that have not yet been reported in the 
observation, but which could become observable in the future 
with better instruments and larger, more complete data bases. 



1.2. Special role of [Fe/H] 

Population synthesis is a powerful tool to investigate such corre- 
lations, because all properties of the (numerical) parent disks, 
and those of the synthetic planet populations emerging from 



them are known. For actual exoplanets this is obviously not the 
case. Observationally, the host stars metallicity only can be de- 
termined and assumed to be a proxy for the disks dust-to-gas 
ratio /d/g- Other fundamental parameters of the disks, such as 
their initial gas mass (which together with /d/g sets the absolute 
amount of solids), and their lifetime, cannot be observationally 
determined. It is therefore no surprise that the observationally 
inferred correlations involve essentially only the metallicity. 

A recently discussed correlation that could partially change 
this is the possible enhanced lithium depletion measured in Sun- 
like stars hosting planets (Israelian et al. I2009I but see Baumann 
et al. (12010b for an opposed view). This is because such a deple- 
tion could result from a positive correlation of disk lifetime and 
likelihood of forming giant planets (Bouvier 2008 ). The influ- 
ence of the disk lifetime Tdisk on the occurrence of giant planets, 
as observed in our models, is discussed in Sect. |8.3| 

In general however, only probability distributions for disk 
masses and lifetimes can be inferred from the observations of 
star forming regions. Our study allows to a certain degree to nar- 
row down these probability distributions for a specific star-planet 
system, as not all types of planets can form in all types of disks. 



1.3. Earlier works 

Several studies have already addressed correlations between disk 
and planetary properties: Ida & Lin d2004bb studied the influ- 
ences of the metallicity as modeled by the dust-to-gas ratio. They 
find that the normalized mass and semimajor axis distributions 
of the potentially observable (giant) synthetic planets are rather 
independent of [Fe/H], while the frequency of giant planets in- 
creases with it, in agreement with observations. Kornet et al. 

(2005) modeled the evolution of the solids from dust size to plan- 
etesimals and, using the final planetesimal surface density, esti- 
mated the giant planet formation capability. Such an approach 
also reproduces the "metallicity effect". Dodson-Robinson et al. 

(2006) derived a fit for the time until gas runaway accretion is 
triggered by a planetary core at a fixed distance of 5.2 AU, as a 
function of the planetesimal surface density. They used this fit 
to estimate the frequency of giant planets as a function of sev- 
eral disk properties. We show our results regarding this subject 
in Sect. [8] comparing with their results. Matsuo et al. (2007) 
used parameterized formation models to determine the disk pa- 
rameters where giant planet formation is possible either by core 
accretion or direct gravitational collapse. We present a similar 
study in Sect. [5] 



1.4. Structure of the paper 

The paper is organized as follows: Section[2]describes the meth- 
ods used to obtain the result concerning the synthetic popula- 
tion presented in Sect. [3] Section[4]studies various aspects of the 
impact of disk properties on the mass of extrasolar planets, in- 
cluding the planetary initial mass function. Section [5] analyzes 
under which disk conditions giant planets can form. Section [6] 
addresses the correlations of disk properties with planetary mi- 
gration and the final semimajor axis distribution. In Sect. [7] we 
study the metallicity of close-in planets. Section [8] shows how 
disk properties determine the fraction of stars with giant planets, 
while Sect, passes the consequences if disk masses and lifetimes 
are correlated. Finally, in Sect. 10 we summarize the results and 
present our conclusions. 
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2. Methods 

As our approach to planet formation and population synthesis 
was described in details in Paper I, we limit here ourselves to a 
short overview. 



2. 1 . General Procedure 

To obtain a synthetic population of planets, we proceed in five 
steps: (1) The probability distributions for the initial conditions 
are derived from observations of protoplanetary disks. (2) A 
large numbers of sets of initial conditions are drawn from these 
probability distributions in a Monte Carlo fashion. (3) The corre- 
sponding final outcomes of the planet formation process (plane- 
tary mass and position) are computed using our planet formation 
model. This results in a population of synthetic planets, most 
of them however undetectable with current observational tech- 
niques. (4) To obtain the subset of the potentially observable 
synthetic planets, we apply an appropriate detection bias. Since 
most planets have been discovered by radial velocity techniques, 
we use a bias based on the velocity amplitude of the host star. (5) 
Correlations between the initial (disk) conditions and the plan- 
etary properties are searched and compared to observations (if 
existing). 

2.2. Initial Conditions - Probability distributions 

We use four Monte Carlo variables to specify the initial con- 
ditions. (1) The dust-to-gas-ratio /d/g which we link to [Fe/H] 
as [Fe/H]=log(/ D/G // D/G , G ), where / D/G;0 is the dust-to-gas ra- 
tio of the solar nebula for which we assume a value of 0.04. 
The choice of this value was discussed in details in Paper I. 
Here we only briefly mention that the factor two to three by 
which 0.04 is larger than the measured photospheric Z of the 
sun (Lodders 2003 ) is a first order representation of the effects 
of dust evolution and drift. These mechanisms lead to an increase 
of the "planetesimal" /d/g compared to the original "dust" /d/g 
in the inner, planet forming parts of the disk by a similar factor 
through the advection of material from the outer disk (Kornet 
et al. 120041) . The probability of occurrence of a given [Fe/H] 
is derived from the metallicity distribution of the FGK stars in 
the CORALIE planet search sample (Udry et al. 120001 1 which 
is representative for the [Fe/H] distribution of the target stars 
in the various major radial velocity search campaigns (Paper I). 
Thus, we assume that the observed stellar metallicity is a good 
indicator of the primordial disk metallicity (Santos et al. 2003 
but see also Pasquini 120071) . (2) The initial gas surface density 
So at «o = 5.2 AU which we link to the initial disk gas mass 

■^disk = 47r£o a o^ 2 ( V a max _ V^min)- The probability distribution 
of circumstellar disk masses is derived from the observations 
of the p Ophiuchi star formating region (Beckwith & Sargent 
1996). (3) The rate of photo-evaporation M w . The distribution of 
M w is constrained by the observed disk age distribution: We have 
adjusted our M w distribution to obtain together with our value of 
the viscosity parameter a a distribution of disk lifetimes 
that is in good agreement with the observed distribution (Haisch 
et al. 120011 ). (4) The initial semimajor axis of the planetary seed 
tfstart- The distribution of the starting positions of planetary em- 
bryos is not constrained by observations and only theoretical ar- 
guments can be used. They indicate (e.g. Mordasini et al. 2009a) 
that runaway bodies should emerge with a uniform distribution 
in log(a stalt ), as already adopted by Ida & Lin (2004a). 

In addition, we assume that these Monte Carlo variables are 
independent variables. Hence, we do not take into account po- 



tential correlations arising during disk formation itself e.g. via 
opacity effects. The fact that we draw Sq and M w independently 
implies, on average, longer disk lifetime Tdi s k for more massive 
disks. The importance of this coupling is discussed in sect. [9] 
We do not consider here the influence of varying initial disk 
radii (Kornet et al. 2005 ). This will be considered in future work, 
taking into account recent observational results (Andrews et al. 
I20091I20T01 ). 

2.3. Planet Formation Model 



We use a slightly modified version (see Mordasini et al. 2009a) 
of the extended core accretion formation model described in de- 
tails in Alibert et al. (2005a). As in Pollack et al. ( 1996), we com- 
pute the evolution of the planetary core and envelope structure, 
but include disk evolution using the a formalism, and planetary 
migration (isothermal type I and type II). We have been able to 
show (Alibert et al. 2005b ) that the model reproduces many ob- 
servational constraints imposed by our own giant planets. 

2.3.1. Relevant model assumptions 

For the disk-planet correlations discussed in this work, a number 
of model assumptions were found to be particularly relevant and 
have directly visible consequences in the correlations. 

The first one concerns the structure of the protoplanetary 
disk, and specifically the location of the iceline. In the nomi- 
nal model, the position of the ice line is an increasing function 
of disk mass due to viscous dissipation (Paper I and III; Min et 
al. 2011). It is found that this has important consequences both 
on the mass of giant planets (Sect. 4.3.1 1 and their formation 
location (Sect. |5.2. 

A second assumption which is directly relevant for the 
masses of giant planets (Sect. |4.3.2| |4.4| ) is that we assume that 
gap formation does no reduce the gas accretion rate of giant 
planets (Lubow et al. 1999). As explained in Paper I, the rea- 
son for this is that Kley & Dirksen (2006) have shown that when 
a giant planet becomes sufficiently massive, the disk-planet sys- 
tem can undergo an eccentric instability. The planet then leaves 
the clean parts of the gap, resulting in a substantial increase of 
the gas accretion rate. This means that the gas accretion rate is 
obtained for low mass planets (M core < 10M ffi ) by solving the 
planetary structure equations while for more massive planets in 
the disk limited, runaway gas accretion phase it is equal to the 
rate at which gas viscously flows towards the star in the disk 
^enve = A^disk = 37rv£ where v is the disk viscosity, and 2 the 
gas surface density. 

A third setting which is relevant for the correlation of migra- 
tion and metallicity (6.1.2 1 is that we assume that the planetes- 



imal accretion rates are independent of the migration rate. The 
M core is calculated in the same way as in Pollack et al. (fT996l 
Possible shepherding effects are not included for the reasons ex- 
plained in Paper I. This means that cores in low solid surface 
density disks still can grow relatively massive by migration. 

Fourth, we assume for type II migration that as soon as the 
planet is more massive than the local disk mass (M p i anet > Sa 2 
where a is the semimajor axis), the migration rate is given as 
-(3v/a) x (Ea 2 /Mpi anet ) p . In the nominal model, we use p — 1 
("fully suppressed" case, Armitage 2007). This is important to 
understand the absence of a strong correlation of [Fe/H] and the 
semimajor axis (Sect. |67T7T||6.3.3| ) 

Note that we simplify the problem significantly by assuming 
that only a single planet can form in a given disk. This would be 
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correct in the limit that proto-planets can form within a single 
system without influencing each other. This is clearly an ideal- 
ization, but in the context of this study this is not necessarily 
a disadvantage. It allows one to see clearly disk-planet correla- 
tions, which otherwise might get partially blurred due to the ran- 
dom character of the (gravitational) interactions between several 
protoplanets. Thommes et al. (2008 ) discuss extensively some of 
the effects induced by the concurrent formation of several plan- 
ets in one disk. 

The population discussed here is essentially obtained, except 
if otherwise mentioned, with identical parameters and Monte 
Carlo distributions as the nominal population presented in Paper 
I and II. In particular, this means that the stellar mass is equal to 
one solar mass M Q , the gas disk viscosity parameter a is 0.007, 
and the type I migration efficiency factor fi is 0.001. 

However, there is one relevant aspect in which the procedure 
used here differs from the one used in paper I and II. We draw 
the initial semimajor axis of the starting seed strictly uniform 
in log(a stalt ) and disregard the additional criteria mentioned in 
Paper I. This eliminates some correlations present already in the 
initial conditions and thus makes it easier to identify the influ- 
ence of disk properties on the planetary properties. 

As a practical unit for the disk gas masses we define, in anal- 
ogy to [Fe/H], a relative logarithmic unit, denoted [Mq/Msn] = 
log (So/£o ; snX where we assume an initial gas surface density 
at 5.2 AU S(),sn = 200 g/cm 2 . This is somewhat more than a 
plain MMSN values - Hayashi's (1981) value would be about 
145 g/cm 2 - but with this choice [Md/Msn] covers, similar as 
[Fe/H], a nearly symmetric range around zero between -0.6 and 
+0.7 for the values of So considered here (50-1000 g/cm 2 , cor- 
responding to disk masses between about 0.004 and 0.09 M Q ). 
As the disk gas mass M D and the initial gas surface density So 
are directly proportional to each other, these two terms are often 
used in an interchangeable way. 



3. Mass-Distance Diagram 

To provide an overview of the synthetic planet population we 
analyze in this work we plot in Fig.[T]the mass-distance diagram 
of the entire population. In the figure, a number of structures al- 
ready discussed in Paper I can be identified again, like the low 
mass "failed cores", the approximately Neptunian mass "hori- 
zontal branch" or the massive "outer group" planets outside a 
few AU. We also see that there are no very massive planets close- 
in. This is a consequence of the fact that planets with a mass 
larger than the local disk mass migrate at a reduced rate because 
their inertia is too large for the finite angular momentum flux in 
the disk. Their absence at small a depends upon the assumed de- 
gree of reduction of the type II migration rate once the planet is 
massive (Syer & Clarke [T9951 

3.1. Paucity of low mass, close-in planets 

There is also a relative paucity of low mass (4 < M < 10M e ) 
close-in (a <; 0.3 AU) planets. Recent observations (e.g. Howard 
et al. I2010t Mayor et al. I201 II) rather indicate a high frequency 
of low mass Super Earth planets close to their parent star. 

Concerning this issue, we first remind that our model shows 
the state of (proto-)planets at the moment when the disk disap- 
pears, and not the final one after billions of years of evolution. 
At this stage, the low mass protoplanets at small distances have 
a mass about equal the local (solid) isolation mass (type I migra- 
tion is strongly reduced in this simulation), which is of order 0. 1 
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Fig. 1. Mass-distance diagram of the synthetic population ana- 
lyzed in this work. The dashed curve shows the feeding limit 
(cf. Paper I): planets reaching this limit have been arbitrarily set 
toO.lAU. 



M @ only. Further growth from these masses to the final masses 
by giant impacts (e.g. Marcus et al. [2010 ) will then only oc- 
cur after the damping influence of the gas disk is gone (Ida & 
Lin |20 10), a phase which is however not included in our model 
at the moment. This effect would populate the depleted region 
from "below". 

Second, the high observed frequency could be indicative of 
quite efficient type I migration, at least in some parts of the disk. 
This would populate the depleted region from "outside". Large 
quantities of low mass, close-in planets are found in some syn- 
thetic populations too, but only for high type I efficiency factors 
(/i ~ 0.1 - 1) as shown in Fig. 11 of paper II. On the other 
hand can the distribution of semimajor axes of giant planets at 
larger distances only be reproduced when the type I migration 
rate (as found from linear theories for isothermal disks) is re- 
duced by a significant factor (jj < 0.01, Paper II; Schlaufman et 
al. 120091) . Here we are using the same low type I efficiency factor 
fi = 0.001 which lead in Paper II to the best reproduction of the 
observed properties of giant planets. Unfortunately, this causes 
in the same time necessarily also the depleted region. 

The answer to this dilemma is probably a significantly more 
complex migration pattern than what can be mimicked with 
global (independent of planet mass and distance) efficiency fac- 
tors like f\ for isothermal type I rates as done here. Recent stud- 
ies of type I migration dropping the often inappropriate assump- 
tion of isothermality (e.g. Paardekooper et al. 120101 Masset & 
Casoli 120101) indeed find complex patterns with rapid in- and 
outward migration. First exploratory planet population synthe- 
ses using such updated type I migration models (Mordasini et 
al. BOllbl ) find that type I migration is directed outward in some 
parts of the disk, and inward in others, which leads to the exis- 
tence of convergence zones (migration traps) which can result in 
a pile-up of many low mass planets in certain parts of the disk 
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Mass [M e ] Mass [M e ] Mass [M e ] 

Fig. 2. Planetary initial mass function at the moment when the gaseous disk disappears. Note that the models doesn't include any 
growth processes like giant impacts after disk dispersal which can significantly modify the mass distribution afterwards, especially 
for low mass planets with M < 10M ffl . Left: PIMF binned according to metallicity. Center: PIMF binned according to the gas disk 
mass. Left: PIMF binned according to the disk lifetime. The meaning of the different lines, and the limiting values, are indicated in 
the plot. 



(cf. Lyra et al. 2010; Sandor et al. 1201 11 1. This important subject 
will be addressed in dedicated work (Dittkrist et al. in prep.). 

A third mechanism occurs when several planets form and 
migrate concurrently. Then, lower mass planets can be pushed 
close to the star after being captured in mean motion resonances 
of a more massive, more rapidly migrating outer planet. This 
will be studied in oncoming simulations allowing the formation 
of many planets per disk (Alibert et al. in prep.) 

4. Mass 

4.1. Planetary initial mass function (PIMF) 

A central outcome of population synthesis is the planetary initial 
mass function PIMF (Paper II and III). Figure [2] shows the PIMF 
of all synthetic planets with a mass larger than 1 M m , binned 
into low, medium and high metallicity (left panel), disk gas mass 
(central panel) and disk lifetime (right panel). The diving val- 
ues are indicated in the plot and are chosen for all three cases 
in a way that the central bin contains about 70% of the plan- 
ets, and the other two about 15% each. We recall that our mod- 
els start with an initial seed mass of 0.6 M @ , so that we cannot 
reliably make predictions about the exact form of the PIMF in 
the < 10M e domain as mentioned in Paper I and II. Note that 
each bin has been normalized individually, so that the absolute 
height of different bins (e.g. low vs. medium metallicity) cannot 
be compared in absolute terms, but taking into account that the 
initial distribution of [Fe/H] follows the distribution observed in 
the solar neighborhood. However, one can directly see the rel- 
ative importance of a given planetary type within a bin (e.g. at 
high metallicity there are more Jovian than Neptunian planets, in 
contrast to the low [Fe/H] case.). One recognizes again a num- 
ber of features discussed in Paper I, like the high mass tail in 
the Super-Jupiter domain, the giant's plateau in the Jovian mass 
regime, the minimum at about 30-40 M @ corresponding to the 
planetary desert (Ida & Lin 2004a), the Neptunian bump and the 
strong rise towards small masses. 



4.1 .1 . PIMF as function of metallicity 

The left panel shows that the metallicity just modifies the im- 
portance of the three basic families of planets visible in the 
PIMF (Jovian planets, Neptunian planets, proto-terrestrial plan- 
ets), without however changing the general shape of the PIMF 
(in contrast to the other two disk properties). For giant planets 
this means that in a high metallicity environment, more giant 
planets form, but their mass distribution is similar. The reason 
for this is related to the fact that a certain [Fe/H] acts as a thresh- 
old for giant planet formation (threshold solid surface density to 
reach a critical core mass, see Sect.|5]l, but is not important in de- 
termining the final total mass, because in the end, gas makes up 
for most of the mass of a giant planet, and not solids. Therefore, 
a high metallicity mainly allows a larger number of high mass 
planets, but not of a higher mass (except for very large masses, 
see sect. 4.3 1. This increase in frequency is of course the under- 



lying reason for the observed metallicity effect, i.e. the increase 
of the detection rate of giant planets with [Fe/H] (see Sec. |8. 1[>. 



Moving down in mass, we see that in the Neptunian mass 
domain, the dependence of planet frequency on metallicity is 
weak. This means that no metallicity effect is predicted in this 
domain. This is in good agreement with recent observations 
(Mayor et al. I201 U . Moving further down in mass to the proto- 
terrestrial domain, the lines are inverted relative to the giant 
planet domain, which means that an inverse metallicity effect 
occurs (low mass planets are more frequent at low [Fe/H] com- 
pared to high [Fe/H]). Concerning this last point, we must take 
into account that in our model, only one embryo can form per 
disk. The situation that in high [Fe/H] disks, both a giant planet 
and a "byproduct" terrestrial planet form is therefore not possi- 
ble. This could artificially strengthen the [Fe/H] -low mass planet 
anti-correlation, so that more secure predictions will be possible 
with models dropping the one-embryo-per-disk simplification 
(Alibert et al. in prep.). Observationally, future high precision 
observations (e.g. with ESPRESSO) will test this prediction. 
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Fig. 3. Planetary mass function binned ac- 
cording to [Fe/H] as in Fig. [2] left panel, 
but now only for the sub-population of po- 
tentially detectable synthetic planets with 
an RV precision of 10 m/s (left panel) 
and 1 m/s (right panel). Green dashed 
lines again represent high, blue solid lines 
medium and red dotted lines low metallic- 
ities, where the threshold metallicities are 
the same as in Fig. [2] 



One further notes that for all three bins, the highest peak in 
the mass function is found for the proto-terrestrial planets, with 
the clearest dominance at low metallicity. For high metallicity, 
the second highest peak occurs for the giant planets, whereas 
for intermediate and even more clearly for low metallicity, the 
Neptunian planets are responsible for the second highest peak. 

4.1 .2. PIMF as function of disk gas mass 

The central panel shows that the disk mass distribution directly 
affects the shape of the PIMF and not only the height of the 
peaks of the distribution as metallicity does. For giant plan- 
ets, a high [Md/Msn] shifts the formation of giant planets from 
lower mass to higher masses (compare the blue and green line). 
For lower disk masses, massive giant planets (> 6Mq^) can- 
not form at all. In this case, there is simply not enough gaseous 
material available. The distribution of intermediate-mass planets 
(30 M® < M < 1 M^), on the other hand, is little affected. We 
conclude that in our model there is a direct correlation between 
the (maximal) mass of giant planets and the disk gas masses. 
This is essentially due to the fact that giant planets accrete most 
of their mass in a regime where their accretion rate is propor- 
tional to the disk mass, see the discussion in Section |4~4| 

In Paper EH, we showed that, in our alpha-disk model, the 
disk mass has to be scaled roughly linearly with M* in order 
to reproduce the observed correlation between the star's accre- 
tion rate and its mass. As shown in this paper, this translated 
into the formation of planets of a higher mass orbiting stars of 
a larger mass, in good agreement with observation (Lovis & 
Mayor 2007 ). Even though we did not vary the stellar mass in the 
present study, the fact that more massive planets form in disks 
taken from the high-mass end of the distribution stems from the 
same reason. Unfortunately, in practice we cannot infer the pri- 
mordial disk mass for any actually detected exoplanet and so 
this correlation is difficult to test observational in contrast to the 
metallicity correlation. The panel further shows that the relative 
importance of Neptune-like planets strongly depends on the disk 
gas mass. Neptune-like planets form particularly easily in disks 
with small primordial disk masses while they do not seem to be 
able to form in massive disks. This again is a consequence of the 
same effect discussed above and already pointed out in Paper III. 

4.1 .3. PIMF as function of disk lifetime 

The third panel at the right of Fig. [2] finally shows the disk life- 
time. The disk lifetime has a twofold influence. In disks with 



longer lifetimes, cores will be able to grow to the critical mass 
and accrete gas in a runaway fashion even for relatively low solid 
surface density. Therefore, the disk lifetime acts as a threshold 
for giant planet formation, similar to the metallicity. The life- 
time of disks however also affects the total mass of the planets 
similar to the disk mass, as, to first order, the planet's final mass 
will be equal to the accretion rate times the duration of the ac- 
cretion phase. This twofold effect is clearly seen in the figure. 
Giant planets formed in long-lived disks are both numerous and 
of a higher mass. Compared to metallicity and disk masses, the 
disk lifetimes has thus a more complex influence on the result- 
ing planet population. The disk lifetime both scales and distorts 
the shape of the PIMF. Similarly to the primordial disk mass, 
it is difficult to deduce the disk lifetime for any give observed 
system. We come back to that in sect. [8] 

4.2. Observable distribution as function of [Fe/H] 

It is interesting to look whether some of the correlations be- 
tween disk properties and the underlying mass function (of all 
planets) discussed in the section above can already be seen in 
the observational data we have today, which represent only a 
small fraction of all existing planets. For this, we plot in Fig. [3] 
the mass histogram binned according to [Fe/H], but now includ- 
ing only planets detectable by a 10-year duration radial velocity 
(RV) survey with a precision of either 10 m/s (left panel) or 1 
m/s precision (right panel). Note that, like as in Paper III, we use 
a very simple velocity amplitude cut-off criterium to determine 
detectability. 

4.2.1. 10 m/s radial velocity precision 

The figure shows that at a precision of 10 m/s, the mass distri- 
bution has a very similar shape for all three metallicity bins, as 
expected from the discussion above and the previous work by 
Ida & Lin (2004b ). However, a closer look indicates that the low 
[Fe/H] bin has a somewhat narrower distribution, and that there 
is a certain systematic difference at the upper mass end (see the 
next section). However, the medium and high metallicity bins to 
which most of the present day known planet population belongs, 
are very similar. It is therefore not surprising that, given the rel- 
atively small number of observed planets compared to the large 
number of synthetic planets used here, no significant correlation 
between the shape of the mass distribution of giant planets and 
the host star metallicity has been noticed thus far. 
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4.2.2. 1 m/s radial velocity precision 

At a precision of 1 m/s on the contrary, the influence of [Fe/H] 
on the PIMF becomes visible as measured by the relative fre- 
quency of Neptunian versus Jovian planets, a trend that has been 
observed (Udry et al. 120061 . We find that the ratio of the number 
of Neptune-mass to Jupiter-mass planets strongly correlates with 
the metallicity. For the low metallicity, the Neptune-mass plan- 
ets are almost as numerous as Jupiter-mass planets while for high 
metallicity disks they are significantly less frequent. This corre- 
sponds well with the observed result (Sousa et al. [2008 ). The 
numerical values of the ratios of synthetic Jupiter-to-Neptune 
number of planets (dividing mass defined at 30 M®) obtained for 
a 10 year search at 1 m/s is in our nominal simulation 4.7, 5.8 
and 7.6 for the low, medium and high [Fe/H] bin respectively. 
These values are difficult directly with Sousa et al. (2008 ) due to 
several reasons (differing [Fe/H] bins, differing primary masses, 
differing survey duration). However, these authors find values for 
the Jupiter-to-Neptune ratio of 1.6-5 for [Fe/H]< -0.15, and 7.5- 
30 for [Fe/H]> 0.15. Note that the ratios even at 1 m/s precision 
are higher than those in the full underlying distribution without 
detection bias. The reason is that even at 1 m/s many Neptune- 
like planets remain undetected. Finally, we point out that due to 
our assumption of a single planet-per-disk simplification, it is 
very likely that our ratio is an upper limit. A hint in the same di- 
rection can also be inferred when comparing our synthetic result 
with the observed one shown in Fig. 10 of Mayor et al. (1201 11 1. 
Otherwise, observed and synthetic result are similar. 

From the two plots we also understand that the decrease of 
the observed mass distribution at 10 m/s towards small masses, 
while being clearly the consequence of the observational detec- 
tion bias against low mass planets, is amplified by the actual true 
decrease of the underlying mass distribution. This degree of de- 
pletion of intermediate mass planets (30 - 100 M e ) is potentially 
an important constraint on the gas accretion rate at the begin- 
ning of runaway gas accretion. In Mordasini et al. (1201 lal l. we 
presented a detailed discussion about various assumptions about 
gas accretion that directly determine this portion of the theoreti- 
cal PIMF. 



4.3. (Maximal) planetary masses and metallicity 

To emphasize the importance of the metallicity (measured in our 
model by the ratio of dust-to-gas in the disk) on the resulting 
high-end of planetary masses, we show in Figure |4]the mass of 
synthetic planets (in units of Jupiter masses) as a function of 
metallicity. While we have indicated in Sect. |4. 1 1 that metallic- 
ity does not change significantly the distribution of the mass for 
the bulk of the population, we see here that the metallicity deter- 
mines the maximum mass a planefjjcan grow to in a given disk, 
in particular for subsolar metallicities. This can be understood 
as follows. To grow to a very high planetary mass, a critical core 
must form before viscous evolution and photoevaporation of the 
disk have had time to significantly deplete the disk mass. Such a 
very early start is not possible in a low metallicity environment 
in which growing large cores takes longer. 

This second order correlation via the core formation 
timescale has been checked and confirmed by artificially reduc- 
ing the starting time of the seed embryos by a factor of two. In 
this test, the maximal mass of planets at [Fe/H] =-0.4 increases 
to about 18 Mq^_ compared to ~ IMq^ in the nominal case. 

1 We here call for simplicity all objects that form in the population 
synthesis planets, even if their mass is larger than the deuterium burning 
limit at about 13 Mq^_. 
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Fig. 4. The mass of synthetic planets (in Jupiter masses Mj) as 
a function of [Fe/H] for the nominal population (small dots). 
There is an absence of very massive planets around low metal- 
licity stars. In the plot, lines indicate approximatively the region 
of high mass and low metallicity where no synthetic giant plan- 
ets are found in other non-nominal populations: with an ice line 
fixed to 2.7 AU (dotted line), or with the effect of gap formation 
on the gas accretion rate modeled as in Veras & Armitage (2004) 
(solid line). The dashed line shows the limit given in Fischer & 
Valenti (125051 



4.3.1 . Effect caused by the iceline position 

A less pronounced absence of massive planets at low metallic- 
ity can also be obtained if we fix the position of the ice line at 
2.7 AU independently of [M D /M$ N ], as expected for an optically 
thin disk irradiated from a 1 M e star (Ida & Lin 2004a). This dif- 
ference is due to an interesting chain of correlations: 

1) We recall that in the nominal case, the position of the ice line 
is an increasing function of disk mass. 

2) At low metallicities, high disk gas masses are needed for giant 
planet formation (Matsuo et al. 120071 Sect. 5.3 1. Thus, the ice 



120071 Sect. _ 

line position in low metallicity disks forming giant planets will 
typically be located at large distances. 

3)Hence, the typical starting position of giant planets-to-be will 
also be large at low [Fe/H] (a sta rt Uk 



5.2 1 



see Sect. 

4)This distant starting position implies a slower growth of the 
core (e.g. Paper I), disadvantageous for successful giant planet 
formation. 

Thus, a positive correlation of the disk gas mass and the iceline 
position renders giant planet formation at low [Fe/H] even more 
difficult, and makes it impossible to populate the upper left cor- 
ner of Fig. [4] with the nominal model. It also means that there 
is a difference between a metal poor-gas rich and a metal rich- 
gas poor environment, even if the solid surface density beyond 
the ice line are in principle the same. A low [Fe/H] cannot in all 
circumstances be compensated a high [Md/MsnL 

In the case with an ice line independent of [M D /Ms N ], this 
complex chain of correlations between [Fe/H], [Md/Msn], ^ice, 
a st art and 4tart is broken, and planets with about 10 can still 
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form at [Fe/H]=-0.5, and that maximal masses become indepen- 
dent of [Fe/H] already for [Fe/H]> -0.3. The upper limit of pos- 
sible masses is indicated for this case by the dotted line in Figure 
[4] It is clear that this chain of correlations depends quite sensi- 
tively on a number of specific model assumptions, making it a 
less robust prediction. 



4.3.2. Very massive planets 

The absence of very massive companions formed by core accre- 
tion at very low metallicities is also seen in the models of Ida 
& Lin d2004bl) and Matsuo et al. (120071 1. However, compared to 
these works, planets of significantly higher mass can form in our 
simulations. This is a direct consequence of the fact that we do 
not limit gas accretion due to gap formation, as mentioned in 
Sect. |2.3.1| If we instead limit in a non-nominal population the 
gas accretion rate by giant planets due to gap formation by us- 
ing the fit of Veras & Armitage (2004), the maximal mass of all 
planets in the population are reduced to ~ IMq^. The absence of 
the most massive objects at low [Fe/H] would remain similar but 
with weaker dependence of the maximal mass on [Fe/H] when 
scaled to smaller absolute masses. The approximative limiting 
envelope for such a population is also shown in Fig.|4]as a solid 
line. This envelope is now similar to the one of Matsuo et al. 
(2007). We thus see that while the tendency towards an absence 
of very massive planets at low [Fe/H] is a general prediction of 
the core accretion theory, the quantitative results depend on spe- 
cific model settings. 

When studying the figure one should keep in mind that very 
massive planets (> lOM^) are in fact very rare outcomes in 
our simulations, in agreement with the observed "brown dwarf 
desert" as discussed in Paper II. Figure|4]may provide the some- 
what misleading impression that these are common objects while 
in fact these objects appear only because the underlying syn- 
thetic population is extremely large (about 200 000 initial con- 
ditions). Even if such a population vastly exceeds the actually 
observed population, these high numbers of planets are required 
to investigate the different correlations some depending on two 
variables (like disk mass and [Fe/H]). 



4.3.3. Comparison with observation 

The absence of very high mass planets (or of a very high to- 
tal mass locked up in planets in multi-planet systems) at low 
metallicities was noted also observationally in Udry et al. (f2002 ), 
Santos et al. ( 120031 and Fischer & Valenti ( |2005T l. The dashed 
line in Fig. [4] is taken from the latter work. It is clear that such 
an absence could also simply be a small number effect: Giant 
planets around low [Fe/H] hosts are rare, and high mass giant 
planets (> SMq^) are rare at all metallicities, which could com- 
bine into the observed paucity. It is interesting to compare our 
results with the observational database especially since the lat- 
ter has been significantly extended since the above mentioned 
papers have been published. 

Such a comparison is presented in Fig. [5] which shows the 
masses of actual and synthetic planets as a function of metallic - 
ity. Observational data was taken from J. Schneiders Extrasolar 
Planet Encyclopedia. Note that this database does not include 
companions larger than about 20 Mq^_ . However, such compan- 
ions are anyway extremely rare, both in the model (only 0.1% 
of the synthetic planets have a mass larger than 20 M^) but 
also observationally. As a result, our observational knowledge of 
this mass range, where several formation mechanism could con- 
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Fig. 5. Observed and synthetic planetary masses as a function 
of metallicity. Black dots are synthetic planets in the nominal 
model, while solid, dashed and dotted lines again indicate the 
limiting envelope of non-nominal models, as in Fig. |4] Red stars 
are observed RV companions around stars with 0.8 < M»/M < 
1 .2. Green squares are observed RV planets around stars with a 
mass larger than 1.2 M Q . For both cases, a point in the middle of 
the symbol indicates that the star is a subgiant or giant. Blue tri- 
angles show planets around HR 8799. Names of relevant objects 
are given in the plot. 



tribute, still suffers from low number statistics (see also Sozzetti 
& Desidera EOlOl l. 

In Fig. [5] red stars symbolize planets detected by the ra- 
dial velocity method around primaries with a mass between 0.8 
and 1.2 M , which are the most relevant cases (in the model 
M* — 1 M Q ). Green squares are companions to stars more mas- 
sive than this, also detected by radial velocity. Stars which have 
evolved off the main sequence (if known) are marked in the fig- 
ure with a dot. Note that stellar evolution could be of relevance 
here, as it might invalidate our underlying assumption that the 
photospheric composition measured today correlates with the 
bulk composition of the disk material at formation. This assump- 
tion could be invalid in case of enhanced heavy element settling, 
which for the sun seems to have already lead to a reduction of 
the present day Z as compared to the primordial Zo by more than 
10% (Lodders|2003). 

Looking at Fig. [5] we see that the bulk of all observation- 
ally detected planets falls into regions of the plot where syn- 
thetic planets are also found. Hence, from this observation we 
could conclude that core accretion can account for almost all 
planets currently known. This conclusion was already reached 
by Matsuo et al. ( 120071 1 . It should be noted that surveys that look 
at very metal poor stars do exist (Santos et al.[2()09). 



4.3.4. Relevant individual objects 

However, eight RV planets orbiting five stars clearly lie outside 
the region where giant planets form orbiting a 1 M star in our 
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nominal model. This group is characterized by a much higher 
fraction of giant stars and of stars more massive than 1.2 M m than 
the full sample. It also has a high multiplicity for giant planets. 
As some of the objects are quite peculiar, it is worth to discuss 
individually some of these objects. 

HD 155358 has two planetary companions of rather low 
mass (< 1 M % ) at 0.6 and 1.2 AU (Cochran et al. |2Q07T >. The 
star has a mass of 0.87 ± 0.07M G and a [Fe/H]= -0.68 + 0.07 . 
HD 155358 is a particular star (Fuhrmann & Bemkopf 2008 ). 
It is a very old, thick disk subgiant with a chemical composi- 
tion very far from scaled solar composition. It is significantly 
enriched in alpha-chain nuclei. [Mg/H] is for example with - 
0.36 much higher than [Fe/H]. What matters for the ability to 
form a sufficiently massive core is the surface density of all con- 
densible elements beyond the ice line which can attribute im- 
portant quantities of matter, in particular the a elements O, Si 
and Mg (Dodson-Robinson et al. 120061 Gonzalez I20091 >. Due 
to the enrichment in a elements, the disk around HD 155358 
was depleted in planetesimals maybe only by a factor 2 below a 

[Fe/H]=0 solar composition disk, equivalent to a [Fe/H] 0.3 

if it were to have a scaled solar composition. At such a value, 
forming the anyway low mass giant planets is certainly possible 
with the core accretion mechanism. 

HD 1 14762 wi th its M sin i = 1 1 M^, a=0.36 AU compan- 
ion (Latham et al. 1989) was discussed in this context already 
by Udry et al. (2002). It has been often considered that the sys- 
tem is seen nearly pole on, so that the companion might in fact 
be a late M dwarf (Cochran et al. |199H Hale |1995l ). An alterna- 
tive hypothesis is based on the finding that HD 1 14762 is chem- 
ically and evolutionary very similar to HD 155358 (Fuhrmann 
& Bernkopf 2008 ) which would mean that its effective surface 
density of planetesimals was much higher than one would infer 
from the low [Fe/H] and scaled solar composition. 

HD 47536 is an old K1II I [Fe/H ] =-0.68 giant which is or- 
bited by one (Setiawan et al. 2003 ), possibly two (Setiawan et 
al. 12008 ) quite massive companions (5 and 7 M^) inside a few 
AU. While Setiawan et al. ([2003 ) originally quoted a possible 
stellar mass of about 1 to 3 M , da Silva et al. (2006) more re- 
cently determined a mass of 0.94 + 0.08M*. With this primary 
mass, the object falls clearly out of the envelope of synthetic 
planets in our model. In particular if its two-planet configuration 
is confirmed, this is an interesting object to study with forma- 
tion models, as one then cannot invoke, due to orbital stability, 
a nearly pole on orientation as a possible explanation. Also for 
the direct collapse model this case is probably not obvious to ex- 
plain due to the small semimajor axis of the planets (about 1 .5 
AU for the inner planet), see e.g. Boley (20091. 

BD+20 2457 is a K2II giant, with an estimated mass of 2.8 ± 
1.5M G , a very low metallicity of [Fe/H] =-1.00 + 0.07 w ith two 
very massive companions in tight orbits (Niedzielski et al. 2009 ). 
Even if the very different primary mass makes a comparison with 
the results here difficult, these companions seem to be far from 
the possible parameter space for core accretion. 

The object around y l Leo A (HIP 50583) was found to have 
a projected mass of about 8.8 Mq^_ (Han et al. [2010). The orbit 
of the companion was however recently astrometrically detected 
by Reffert & Quirrenbach (120111 ) in the Hipparcos data. Their 
results indicate a clearly larger actual mass of about 66 Mq^ (with 
substantial incertitude). This moves the companion out of the 
relevant mass domain. 

The plot also includes HR 8799 with a measured [Fe/H]=- 
0.47 (which might however not reflect the initial metallicity, 
Marois et al.|2p08 ) and planets detected by direct imaging i.e. at 
large orbital distances (Marois et al. l2008ll20T0l ). HR 8799 also 



probably has a mass larger than 1.2 M Q (1.5 ± O.3M ), so the 
synthetic population does not apply directly, but it is neverthe- 
less interesting to note that these planets come to lie in a region 
in the [Fe/H] -mass plane which is at or close to the limit where 
giant planets can originate from core accretion. This is a finding 
independent from the fact that the large semimajor axes of these 
planets makes core accretion as the formation mechanism diffi- 
cult (Dodson-Robinson et al. 2009), if no additional mechanism 
causing outward displacement like scattering is acting. 

4.3.5. Summary concerning the [Fe/H]-M correlation 

In summary we see that there are only extremely few examples 
of bona fide, massive companions (5 < M < 20 M^) orbiting 
solar-like, main sequence stars at small orbital distances which 
do not fall into the [Fe/H]-M parameter space covered by our 
implementation of the core accretion model. 

The solid line in Fig. [5] shows also the limiting envelope of 
synthetic planetary masses derived in the case where the gas ac- 
cretion rate due to gap formation is limited by using the simple 
one parameter fit of Veras & Armitage ( 120041 ) to hydrodynam- 
ical simulations of Lubow et al. (1999|. It is clear that such a 
simple fit can only be a rough approximation of the real effect 
and that e.g. disk viscosity also has an influence on the degree 
of quenching of gas accretion (Lissauer et al. 120091) . But the sig- 
nificant number of planets lying above the solid line indicates 
that mechanisms like the eccentric instability (Kley & Dirksen 
2006 ) that allow growth beyond the gap barrier seem to play an 
important role in nature. 

We finally note that the most massive synthetic planets could 
also, at least partially, be an artifact of the one embryo per 
disk approximation. In a disk where several protoplanets con- 
currently grow, planets compete for gas and ejection can even 
remove some planets (Thommes et al. 2008). In this sense, look- 
ing at the total mass in the system should be a more adequate 
quantity than individual planet masses. 

4.4. Planetary masses and disk mass 

Figure [6] shows planetary masses as a function of the initial disk 
gas mass, again in units of Jupiter masses to focus on the massive 
planets. Here we find a different behavior than for [Fe/H]: As ex- 
pected from the discussion of the influence of [M D /Ms N ] on the 
PIMF, we see that there is a positive correlation of disk mass 
and planetary mass over the full domain of disk masses consid- 
ered here. This is illustrated by the dotted line corresponding to a 
mass that is equal half the initial disk mass. This shows that there 
is, for the domain of [Mq/Msn] Z -0.2 an approximatively lin- 
ear correlation between the maximal planet mass and the disk 
mass. The reason for this is that giant planets accrete in our 
model (see Mordasini et al. :2009a) most of their mass in the disk 
limited accretion regime where the planetary gas accretion rate 
is assumed to be the same as in the disk (M env = Mdisk = 37rvZ), 
which is itself proportional to the initial disk mass for a given 
distance and moment in time, at least if the disk evolution is 
close to self-similarity solutions (Hartmann et al. 1998). 

The different influence of [Fe/H] and disk gas masses on the 
mass of giant planets can be quantified: For planets with a mass 
larger than 100 M®, one finds that the median planetary mass 
grows for an increase of [Fe/H] from -0.45 to 0.45 only by about 
a factor 1 .5 (from 5 12 to 747 M @ ). For an increase over the same 
[Md/Msn] domain in contrast, the median giant planet mass in- 
creases from 182 to 814 M e i.e. by about a factor 4.5. 
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Fig. 6. The mass of synthetic planets (in Jupiter masses) as a 
function of the initial disk mass for the nominal population 
(black dots). The dotted line indicates a mass half as large as 
the initial gas disk mass, showing the linear correlation of disk 
mass and maximal planet mass for most of the domain. The solid 
line shows the upper limiting envelope in a non-nominal popula- 
tion where the limiting effect of gap formation on the gas accre- 
tion rate is taken into account, using the fit of Veras & Armitage 
(2004). 



For [Md/Msn] S -0.2, the maximal masses decrease faster 
then linearly with decreasing [M d /Msn] and are thus smaller 
than indicated by the dotted line. This is because at such low 
disk masses, the low resulting solid surface £s = /d/g^o den- 
sities becomes important in a second order effect, so that the 
final planetary mass is no more just determined by the amount 
of gas that can be accreted, but also by the time needed to form 
a supercritical core. This is an effect that is analogous to the de- 
pendence of the maximal mass seen at low [Fe/H]. Note that the 
exact value of the highest efficiency of converting disk gas into 
planetary envelope material (here about 0.5) depends on the ex- 
act treatment of the back-reaction of the planet's accretion on 
the disk, but is found to lie in a domain of 0.3 to 0.5. The lin- 
ear correlation between disk mass and maximal planetary mass 
is however not affected by different treatments. This indicates 
that a positive correlation of disk mass and giant planet mass is 
a stable prediction of the core accretion theory, while the quanti- 
tative degree is model dependent. For the large majority of giant 
planets (M > 1 Mo_) the efficiency of converting disk gas into 
planetary material is lower and to order of magnitude 0.1, but 
with a very wide spread of possible values, depending on e.g. 
the photoevaporation rate, the metallicity or the disk gas mass 
itself. 

The solid line in Figure [6] indicates again the approximative 
upper envelope obtained for the non-nominal population using 
gas accretion rates limited because of gap formation according 
to the Veras & Armitage (2004) fit. In this case, the planetary gas 
accretion rate is still proportional to Mdisk but decreases expo- 
nentially as exp(-M p i anet /1.5 M^) with planet mass until a floor 



value of M env = 0.04Mdisk is hit. The consequence is as shown 
by the figure that the clear correlation of disk mass and planet 
mass is now broken, as expected from the functional form of the 
equation: One finds that the maximal mass now depends only 
weakly on the disk mass, approximatively as (Md/Msn) 1 ^ 4 for 
Mq > Msn- A weak correlation remains because of the floor ac- 
cretion rate. At lower disk masses, the correlation is somewhat 
stronger, for the same reason as for the nominal population. 

The fact that for the nominal population, maximal masses are 
linearly proportional to disk masses is interesting as one there- 
fore expects to see an imprint of the disk gas mass distribution 
on the planetary mass distribution. This is reminiscent of the sit- 
uation for stars, where the protostellar core mass function CMF 
has a remarkably similar functional form as the stellar IMF (e.g. 
McKee & Ostriker 2007 ). This question will be addressed fur- 
ther in forthcoming work. 

5. Disk conditions leading to giant planet formation 

5.1. Minimal solid surface density 

With the population synthesis calculations, one can study a pos- 
teriori which combinations of disk properties allow the forma- 
tion of giant planets, and in particular which are the most ex- 
treme ones. Comparable studies have been done in Kornet et al. 
(120061 ). Dodson -Robin son et al. ([20061 ), Ikoma et al. d2000l or 
Thommes et al. d2008l . 

Figure [7] shows the relative planetesimal surface density 
Ss/2s,sn at 5.2 AU, where 2s = /d/g^o is the planetesimal 
surface density in a given disk and Xs.sn = /d/g,q£o,sn = 8 
g/cm 2 is the planetesimal surface density in our solar nebula 
reference disk, which lead to the formation of a giant planet 
(M > 300M ffl ), as a function of the starting position of the seed 
embryo. The core accretion rate M core °c Xs (e.g. Alibert et al. 
2005a), therefore this quantity is within the core accretion for- 
mation paradigm a central control parameter. Blue filled circles 
represents synthetic planets of the nominal population in all four 
panels. A value of e.g. £s/£s,sn = ^ tnus corres ponds to a disk 
with roughly three times more solids than the amount we as- 
sume for the solar nebula. In all disks, the solid surface density 
scales as a~ x with a jump of a factor 4 at the ice line, therefore 
this value corresponds for a starting position of the embryo of 
e-g- Qstart = 10 AU (which is outside the ice line for all So con- 
sidered here) to a local planetesimal surface density at a start of 
about 3 x 8 x (10/5.2r L5 = 9 g/cm 2 . 

5.1 .1 . Sweet spot for giant planet formation 

The plot shows that the sweet spot for giant planet formation 
(i.e. the lowest solid surface densities which allow that) occurs 
at about 5 AU, at a value of about 0.75, corresponding to 6 
g/cm 2 . This value is similar to the result of Dodson-Robinson 
et al. (I2006I ) who find that a solid surface density of about 6.5 
g/cm 2 (also at 5.2 AU) is needed to bring an embryo to run- 
away gas accretion in 7 Myr, which corresponds to the longest 
living disks in our population (Paper I). At a starting position 
of 1 AU (which is inside the ice line), a £s/2s,sn of at least 10 
is needed, corresponding to a local planetesimal surface density 
of about 0.25 x 10 x 8 x (1/5. 2)' 1 5 ~ 24 g/cm 2 , which is in 
good agreement with Kornet et al. (2006). We see that both at 
smaller and at larger distances, more massive disks of planetesi- 
mals are needed. This is in agreement with the works mentioned 
earlier. The abrupt increase by about a factor 4 at about 3-4 AU 
simply corresponds to the increase necessary to compensate the 
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Fig. 7. Relative planetesimal surface density Zs/Ss.sn = /d/g^o/(/d/g,o2o,sn) at 5.2 AU leading to the formation of a planet with a 
final mass M > 300 M & , as function of the starting position of the embryo a s tart- A value of Es/Es,sn = 1 corresponds to an initial 
planetesimal surface density of 8 g/cm 2 at «o = 5.2 AU, while a value of for example 3 corresponds at e.g. 10 AU to 9 g/cm 2 . 
The top left panel shows the nominal population only (blue dots). The three other panels show non-nominal population as open red 
squares which overlay the nominal population that is also shown for comparison, still with blue dots. Top right: Population with 
faster type I migration (/i = 0. 1). Bottom left: Population where f start is reduced by 50%. Bottom right: Population where the opacity 
in the planetary envelope is 2% of the interstellar value. 



decrease of the planetesimal surface density inside of the ice line 
by an identical factor due to the sublimation of ices. The upper 
limit of Zs/£s,sn m tne pl ot corresponds to disks where both 
fu/o and So come from the upper end of their distributions. 

5.1.2. Isolation mass and timescale effect 

The reasons for the remaining increase of the necessary £s with 
both increasing and decreasing a start i.e. the existence of the 
sweet spot has been discussed in earlier work (e.g. Kornet et 
al. I2006I Thommes et al. l2008l . Therefore we here only briefly 
illustrate this by the following two non-nominal populations 
which are also plotted in Fig. [7] 

The red empty squares in the top right panel correspond to a 
population with a faster type I migration rate (fi =0.1 instead 



of 0.001). In this case, the minimal necessary planetesimal disk 
mass is reduced at small distances. This shows that at small dis- 
tances (where the core accretion timescales are short) the small 
isolation masses (and the associated long Kelvin-Helmholtz 
timescales for gas accretion) limit giant planet growth. Isolation 
masses decrease with decreasing semimajor axis which must be 
compensated by an increasing Ss. Faster type I migration par- 
tially neutralizes this need, as it allows growth beyond the isola- 
tion mass and up to the critical mass by increasing the embryo's 
feeding zone (Alibert et al. 120041 1. 

The empty red squares in the bottom left panel show the re- 
sult for a population for which we have arbitrarily reduced f start 
by a factor 2. In this case at large distances a difference is seen, 
while in the other parts of the plot, the points lie on top of each 
other. At large distances, high amounts of solids are available in 
the embryo's feeding zone, but the core growth timescale is very 
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long. To compensate that, i.e. to still build up a critical core dur- 
ing the disk lifetime, Es must also increase towards the exterior. 
When we (artificially) speed up the accretion rate (or, equiva- 
lently, reduce f st artX this requirement is relaxed. 

5.1 .3. Situation of the solar nebula 

The figure also shows that the solar nebula is not far above the 
overall minimal necessary Xs for giant planet formation. This 
conclusion is also in agreement with Thommes et al. (2008]). We 
find that the overall minimal value is about £s/Xs.sn = 0.75, 
i.e. 2s ~ 6 g/cm 2 . This would correspond to an isolation mass 
at 5.2 AU of about 5.3 M ffi . Studying the total mass of accreted 
planetesimals (which can be in the core, or dissolved in the en- 
velope) we however find that the smallest amount of heavy el- 
ements in planets more massive than 300 M e is clearly larger, 
namely about 23 M e . This is a consequence of migration, and of 
planetesimal accretion also after runaway gas accretion is trig- 
gered. Internal structure modeling of observed transiting extra- 
solar planets by Miller & Fortney (201 1 ) points to a similar min- 
imal heavy element content of giant planets. 

5.1 .4. Impact of the grain opacity 

The bottom right panel of Figure [7] finally shows a population 
calculated with an opacity of 0.02 times the nominal (interstel- 
lar) value in the envelope (e.g. Pollack et al. 1996). Here we find 
that the minimal necessary value is, as expected, lower, namely 
about 4.7 g/cm 2 . This reduction is rather small compared to the 
large effect of the opacity in formation calculations without mi- 
gration (Pollack et al. |1996i Movshovitz et al. l2.Q101 >. The reason 
for the rather small influence is that thanks to migration (even 
with the strongly reduced type I migration rate), cores never 
get completely cut from a supply of fresh planetesimal to ac- 
crete. Therefore, the need for a rapidly growing gaseous enve- 
lope (where the opacity matters) in order to expand the solid 
feeding zone is not as important as in the case of a strict in situ 
formation. Such a growth mode that is limited only by the accre- 
tion of solids, and not the Kelvin-Helmholtz timescale of the en- 
velope, is somewhat related to the situation for Saturn discussed 
in Dodson-Robinson et al. (1200 8). In agreement with this work 
we see that the reduction of the opacity is important for smaller 
distances < 7 AU. 



5.2. Starting position 

We have seen in the last section that the location of the ice line 
is important for the conditions necessary for giant planet forma- 
tion, and influences in this way the birth place of giant planets. 
In Fig. [8] we directly plot the starting position of embryos a sta rt 
which later become planets larger than 300 M & , relative to the 
position of the ice line a la , in the corresponding disk, as a func- 
tion of [Fe/H]. Additionally, the disk mass is color coded. Note 
that in our simulation, the location of the ice line is independent 
of [Fe/H] (but dependent on So). 

The figure shows that at low metallicities [Fe/H]< -0.3, 
planets can only start to form in a small region with a width 
of 2-3 AU just outside the ice line. In this low [Fe/H] disks, 
only in this region massive cores can form quickly enough. 
With increasing [Fe/H], the zone expands to larger radii, and 
for [Fe/H]> -0.15, giant planets can also form inside the ice 
line. This corresponds, with the highest So of the distribution 
(which is necessary, see next section) to a solid surface density 
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Fig. 8. Starting position relative to the ice line of embryos grow- 
ing eventually to planets with a final mass M > 300M ffi as a 
function of [Fe/H]. Symbols show the disk gas mass: Red filled 
squares are low [Md/Msn] < -0.1. Black triangles are interme- 
diate [Md/Msn] (-0.1 to 0.4), and green open circles finally are 
those disk with a large mass [Md/Msn] > 0.4. 



about 3.5 times as high as Ss ; sn ( as a ^ so visible in Fig. [7 1. These 
results are qualitatively similar to those of Ida & Lin (2004a), 
except for a somewhat lower numerical value (3.5 instead of 5). 
The absence of red squares indicating disks with a low gas mass 
at negative values of a start - ai ce (i.e. inside the ice line) shows 
that the disk masses must always be rather large for the forma- 
tion of giant planets inside the ice line, even if [Fe/H] is large. 
At higher [Fe/H]> 0.3, giant planets can form at almost all semi- 
major axes (about 1 AU to 20 AU) if concurrently [Md/Msn] is 
high as indicated by the green open circles, but the preferred for- 
mation location is still beyond the ice line. Finally, at the highest 
[Fe/H]> 0.4, the location of the ice line becomes less important. 

One therefore sees from the figure, that while at low metal- 
licities, pathways to giant planet formation are very restricted 
and in particular only possible in a small part of the disk, this 
is not the case for higher metallicities, where giant planets can 
form all over the disk. 



5.2.1. Anti-correlation of [Fe/H] and a start 

There is another important chain of correlations regarding the 
starting positions, which cannot be directly seen in Fig. [8] but 
which was already mentioned in sect. |4. 3.1] For giant planet for- 
mation at low [Fe/H], high [M d /Msn] are needed for compensa- 
tion (see below). At high [Fe/H] in contrast, low [Md/Msn] are 
sufficient. As the distribution of disk gas masses peaks at values 
that are close to the lower limit of values allowing giant planet 
formation, the typical [Md/Msn] which leads at a high metallic- 
ity to the formation of a giant planet is a rather low [Md/Msn] ■ 
As seen in the figure, also at higher [Fe/H], the ice line remains 
to be the typical formation location. This means that if low disk 
masses correspond to low values of (as we assume in the 
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Fig. 9. Plane of initial conditions [Fe/H]-[M D /M S N] leading to 
the formation of synthetic planets larger than 300 M @ . Different 
symbols indicate disk lifetimes: Red filled squares are low 
Tdisk < 2 Myr. Black triangles are intermediate Tdisk (2 to 4 Myr), 
and green open circles are disks with a long lifetime of 4 to 
7 Myrs. The dotted line corresponds to an initial solid surface 
density of 6 g/cm 2 at 5.2 AU. Between the dashed and the dotted 
line, no short living disk can lead to giant planet formation. 



nominal case), there is a negative correlation of [Fe/H] and the 
typical a start of giant planets-to-be. So in terms of absolute val- 
ues of a start (not relative to the ice line), even if giant planets also 
come in higher [Fe/H] environments typically from beyond the 
ice line, their starting positions are smaller in absolute terms, as 
the typical smaller. The magnitude of this effect can be 

quite significant, and corresponds to a typical difference over the 
1 dex interval [Fe/H] covers of about 3 to 4 AU. If a; ce is in con- 
trast independent of So (which could be possible due to a dead 
zone which we do not include here, or if the disk is optically 
thin), this effect does not exist. It is therefore a finding sensitive 
to model assumptions. 



5.3. Compensation effects between disk mass and [Fe/H] 

Figure [9] illustrates several compensation effects in the metal- 
licity versus gas mass plane, with disk lifetimes coded with 
different symbols. Each point in the graph corresponds to an 
initial condition which allows the formation of a giant planet 
(M > 300M e ). Therefore, this plot shows three dimensions of 
the four dimensional parameter space of initial conditions (astart 
is here projected into the plane). If we would plot all initial con- 
ditions of the population, the full area of the figure would be 
filled with points, distributed as a two dimensional gaussian dis- 
tribution with means just below zero dex both for [Fe/H] (-0.02) 
and [M D /M SN ] (-0.04). 

The figure shows that only initial conditions lying in the up- 
per right corner allow the formation of giant planets. The diago- 
nal line corresponds to a solid surface density of Es = /d/gEo =6 
g/cm 2 , the minimal necessary value identified in Fig.pl The fact 



that the points nearly fill the area above the line show that in this 
region, any combination of the dust-to-gas ratio and the disk gas 
mass allows giant planet formation, provided that the product of 
these two quantities gives a value above the threshold value. This 
means, a low [Md/Msn] can be compensated by a high [Fe/H], 
and vice versa. A certain deviation of this purely multiplicative 
behavior can be seen at very low [Fe/H]. There, the mechanism 
plays that the very high necessary [Md/Msn] come along with 
a so large a lce and thus typically large starting position a sta rt and 
in turn so long core formation timescale, that the necessary Zs 
increases stronger than in the rest of the parameter space, as dis- 
cussed earlier in Sect. 14.31 



5.3.1. Comparison with observations 

It is interesting to compare this result with the works of Greaves 
et al. d20071 > and Wyatt et al. d20071 >. These authors derived, 
purely based on observational data (metallicities of giant planet 
hosts and dust disk mass distribution measured by submillime- 
ter observations) that a minimal mass of solid elements of 
0.5 + 0.1 M^ in a disk is necessary form giant planets. In the 
model, which is in contrast built on first theoretical principles a 
minimal necessary Es of 6 g/cm 2 corresponds to an initial mass 
of planetesimals in the computational disk of about 0.4 to 0.5 
Mq^_ (the spread comes from the position of the ice line which is 
not fixed by the product /d/gEo)- This is in excellent agreement 
with the observed value. In the model, there is no threshold solid 
disk mass of any kind built in. In our models, and in contrast 
to Ida & Lin d2004al I2004bl [20051 120081 . not even a minimal 
core mass for the start of runaway gas accretion is built into the 
model. This threshold mass is not specified, but obtained by the 
solution of the planetary structure equations. This agreement is 
therefore a good indication that our core accretion model catches 
important mechanisms occurring in giant planet formation. 



5.3.2. Long necessary T disk at low solid surface densities 

The colors and shapes of the symbols in Fig. [9] show that there 
is another correlation between the initial conditions, explicitly 
not discussed by Greaves et al. d20071 l: The absence of red filled 
squares (showing short Tdisk) roughly between the dashed and 
the dotted line shows that near the solid surface density thresh- 
old, only long living disks can produce giant planets. This is im- 
mediately understandable within the core accretion model: Near 
the threshold value of Es, core formation is slower than at higher 
solid surface densities, which can be compensated by long disk 
lifetimes to still be able to form giant planets. Note that the distri- 
butions of disk lifetimes of our synthetic disks is by construction 
in agreement with the observed distribution (Paper I). It is inter- 
esting to note further that if the solar nebula indeed had a [Fe/H] 
and [Md/Msn] close to zero dex, we can deduce from Fig. [9] that 
the solar nebula had a rather long disk lifetime of probably more 
than roughly 3-4 Myr, necessary to from Jupiter in such condi- 
tions. This in turn is in good agreement with constraints coming 
from completely different measurement, namely that ages of CB 
chondrules and CV CAI which indicate a minimum lifetime of 
the solar nebula of 4.5 ± 0.8 Myr (Scott 120061 . 



6. Semimajor axis 

In the last section, we mainly studied the effects of disk prop- 
erties on planetary mass. Another important observable quantity 
is the final semimajor axis. As explained in Paper II, one should 
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Fig. 10. Extent of migration Aa — Cl — ^ start ^ 

function of the 

disk gas mass [Mp /Msn] for synthetic planets eventually becom- 
ing giant planets (M > 300 M e ). Symbols indicate disk metal- 
licities: Red filled squares are low [Fe/H]< -0.2. Black triangles 
are intermediate [Fe/H] (-0.2 to 0.2 dex), and green open circles 
are disks with a high [Fe/H]> 0.2. A number of lines are plotted 
to show relevant limiting regimes (see text). 



keep in mind that the results concerning the semimajor axis are 
likely less firm than those for the mass. The reason is first that the 
type I migration model we use here is still based on the isother- 
mal Tanaka et al. (2002) migration rates, which were the only 
existing analytical description of type I at the time the model 
was written. Second we do not include any dynamical effects 
between multiple embryos. Nevertheless, we can gain insights 
of the role of several disk properties on planetary orbits in the 
idealized case of a single (massive) planet growing per disk in 
the limit of nearly absent type I migration. 

6.1. Extent of migration as function of disk mass and [Fe/H] 



Figure 10 shows as a function of disk gas mass (or equivalently, 



of the gas surface density) how far giant planets-to-be (final plan- 
etary mass larger than 300 M e ) migrate, i.e. the difference Aa 
between the final semimajor axis of the planet a and the initial 
semimajor axis a stal1 where its seed embryo is put into the disk. 
Both type I and type II migration rates for the planet dominated 
regime, which is the most relevant one (Paper I), are proportional 
to the gas surface density. Note that by construction, only inward 
migration is possible in the current model (Alibert et al. l2005at : 
For type I migration, the equations of Tanaka et al. (2002]! in 
practice always yield negative torques, and for type II we implic- 
itly assume to be well inside the radius of maximum viscous cou- 
ple, which is justified in almost all circumstances (Paper I). The 
colors and symbols additionally show the corresponding metal- 
licity of the disk. A quite complex pattern is seen, which is the 
consequence of several effects: 

One expects that at a high disk gas mass, the largest extent 
of migration happens. This can indeed be the case, as indicated 



by the dotted line. At the highest [Md/Msn], the extent of mi- 
gration can be very large, up to -8 AU as seen in the lower right 
corner of the graph. But it is more complex: when the disk mass 
is high, the extent of migration can also be minimal as shown by 
the green points in the upper right corner, where only about 0.5 
AU of radial displacement occurs. The difference comes from 
the metallicity, as shown by the symbols. In the bottom right cor- 
ner, there are a many red squares, indicating low [Fe/H]< -0.2, 
which do not exist in the top right corner. (Ignore for the moment 
the few green open circles found in the outermost bottom right 
comer between the dashed and the dotted line. The very large 

displacement 8 AU there are special cases addressed further 

down). 



6.1 .1 . Braking effect 

The reason is that metallicity first acts through a "braking ef- 
fect": When both the gas mass and the metallicity is high, the 
planetesimals surface density is high. Therefore, even if the 
torques would be strong, planets get quickly so massive (un- 
dergo quickly runaway gas accretion) that the disk soon cannot 
force them to migrate on the viscous timescale any more, as the 
planet mass is larger than the local disk mass (Lin & Papaloizou 
1986). Rather, these planets quickly get into the slower planet 
dominated type II migration ( Armitage 2007 ) where the migra- 
tion rate is inversely proportional to the planet mass. Therefore, 
if both the gas mass and the metallicity is high, planets do not 
migrate much if the planetary embryo starts at a rather small 
semimajor axis (which is possible in metal rich disk as we have 
seen above) because there the local disk mass is small (Paper I). 
We see thus that via the planet mass, there exists a negative cor- 
relation between metallicity and migration, even if no influence 
of the metallicity on the migration rate is directly included in the 
model. 



6.1.2. Collection effect 

Second, metallicity acts through a "collection effect": When the 
gas mass is high, but [Fe/H] is low, then the extent of migration 
is large (red squares in the bottom right corner). Under such disk 
conditions, the planetesimal surface density is only intermedi- 
ate, but the torques are strong. Cores cannot become supercriti- 
cal for gas runaway accretion in situ. In order to become a giant 
planet, they must first migrate inwards through the disk, collect- 
ing the planetesimals they come across, until they reach a mass 
where gas runaway accretion sets in (see the formation tracks 
in Paper I). Then the planet slows down due to the mentioned 
"braking effect". In other words, a lower planetesimal surface 
density and a high gas mass make that planets migrate strongly, 
because they stay longer in the faster disk dominated type II mi- 
gration. Therefore, the extend of migration is large, of order -4 
to -6 AU. The collection effect makes that the minimal extent of 
migration (upper limit of the envelope indicated by the dashed- 
dotted line) increases with decreasing disk mass. Note that the 
"braking" and the "collection" effect mean that accretion and 
migration are two interdependent mechanisms that should not 
be treated independently. 

Third, one can see in the graph again the "compensation ef- 
fect" discussed above in Sect. 5.3 i.e. the fact that at low gas 
masses, high [Fe/H] are necessary to form a giant planet. This 
is illustrated by the absence of red squares at low [M D /Ms N ]. 
Instead, if the gas surface density is low, [Fe/H] must be high 
(only green and black symbols). At these rather low disk masses, 
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migration is of intermediate importance, as the torques are not 
very strong, but at the same time the cores have to migrate over 
a certain distance until they have collected enough solids to get 
into gas runaway accretion, and slow down. Therefore, the ex- 
tend of migration is intermediate at low [Md/Msn], about -3 AU. 

Fourth, we finally also notice a consequence of the 
"timescale effect" mentioned before in Sect. 15.1.21 The fact that 
in the outermost lower right corner (between the dashed and the 
dotted line) the symbols are all green circles only is the con- 
sequence of the following: to be able to migrate over such a 
large distance, planets must already start at large initial semi- 
major axes and nevertheless manage to grow supercritical on a 
timescale shorter than the disk lifetime. This is however only 
possible for special cases when the solid surface density is very 
high, otherwise the formation timescale for a massive core is too 
large. Therefore, the planets that migrate most are again found in 
disks with a large mass and [Fe/H], and with a large a stm . They 
don't stop due to the "braking effect", as the local disk mass 
which is given approximately as Y,a 2 is large at large distances. 

6.1.3. Dependence on model assumptions 

We have studied how these results depend on model assump- 
tions, by looking at several non-nominal populations. For a pop- 
ulation where the ice line is fixed at 2.7 AU independently of 
[Md/Msn], the result remains quite similar. 

In a non-nominal population where type II migration is only 
partially suppressed in the planet dominated regime, i.e. where 
the reduction of the type II migration rate due to the planet's in- 
ertia is assumed to scale proportional to (Za 2 /M p i- dnet )^ 2 instead 
of linearly (cf. Alexander & Armitage 2009 ), the triangular en- 
velope shape of the points is approximately retained, but shifted 
to larger extents of migration. For example, the extent of migra- 
tion at the lowest [M d /Msn] increases from -3 to about -5 AU, 
and the maximal overall distance is about -15 instead of -8 AU. 

Using a type I migration efficiency factor of fi — 0.1 has a 
remarkably weak influence on the dependence of Aa on the disk 
properties, except from shifting the whole envelope of points to 
larger values by about -0.5 AU. 

The plot also indicates that for solar nebula like conditions 
([Fe/H] and [M D /M SN ] » 0), a Aa of about -4 AU is found, 
compatible with our earlier results for the formation of Jupiter 
(Alibertet al. l2005bb . 

6.2. Extent of migration as function of Tdisk 



In figure 1 1 we have plotted the extent of migration as a function 



of the life time of the protoplanetary disk. The symbols indicate 
the relative amount of solids in the disk, £s/^s,sn> cf- Fig- 
The boundaries for the three bins were chosen in a way that each 
of the three bins contains about the same number of planets. As 
expected (Alexander & Armitage 2009) disks with a long life- 
time drive migration over larger radial extents. The plot shows 
that the bulk of the planets roughly lies on a line (indicated in 
the plot) running from Aa m -2 AU for the disk with the short- 
est life time which still can form a Jupiter mass planet (about 
0.5-1 Myr) to a Aa a -5.5 AU for the disks with the longest 
Tdisk °f about 7 Myrs. Such a relationship is expected, because 
there is first a characteristic locus from where giant planets come 
(sect. |5.2| i, and second the longer the disk lifetime, the longer 
the torques can act to move the planet. Cases falling on the line 
mostly come from the intermediate bin of 2 < 2s/^s,sn ^ 4, 
which corresponds to typical values of the solid surface density. 
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Fig. 11. Extent of migration Aa — a - a start as a function of the 
disk life time Tdisk for synthetic planets eventually ending up as 
giant planets (M > 300M ffi ). Symbols indicate relative surface 
densities of planetesimals 2s/^s,sn: Red filled squares are low 
£s/2s,sn < 2- Black triangles are intermediate Ss/^s,sn (2 to 4), 
and green open circles are disks with a high Es /Es.sn > 4. The 
solid line indicates the very roughly linear behavior of Aa as a 
function of Tdisk- There are however large systematic departures 
from this simple dependence. 



There is however a significant spread around the diagonal 
line. There are, for example at each Tdisk points which lie above 
the plotted line, i.e. where only very little migration occurs. The 
symbols in this area (green circles and black triangles) indicate 
that these are disks with a high solid surface density. One then 
finds that such cases all correspond to planets starting inside the 
ice line, which is possible for high Es/^s,sn ^ 3.5, see Fig. [7] 
As explained in the previous section, these planets don't migrate 
much, because of the breaking effect that sets in soon. 

On the other hand, there are also cases where clearly more 
migration happens than for the bulk of the planets, i.e. which lie 
clearly below the line. The symbols indicate that this is a mixture 
of both low (red squares) and high (green circles) solid surface 
densities. The former class corresponds to planets that start not 
far outside the ice line, but then migrate all the way down close to 
the star (0.1-0.7 AU), because of the "collection effect", whereas 
the later correspond to planets starting at very large semimajor 
axes, well outside the ice line, which is possible in such metal 
rich disks (cf. Fig. [8}, and migrate then a lot because of the brak- 
ing effect setting in late at such distances, but still never end up 
inside ~ 3 AU. 



6.3. Final semimajor axis distribution of giant planets 

As mentioned in the introduction, there is no clear imprint of 
[Fe/H] onto the semimajor axis distribution of giant planets in 
the observational data, see e.g. Valenti & Fischer, (2008]) or Udry 
& Santos J2007I). I t was studied repeatedly (e.g. Sozetti [2004; 
Santos et al. 12006) whether stars hosting a Hot Jupiter are par- 
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Fig. 12. Distribution of final semimajor axes of synthetic planets larger than 300 M®, as a function of metallicity. In both panels, the 
red dotted lines are for [Fe/H]< -0.2, the blue solid lines for -0.2 <[Fe/H]< 0.2, and green dashed lines are for [Fe/H]> 0.2. The 
left panel is the nominal population, while the right panel is a population where a\ ce = 2.7 AU fixed for all initial conditions. Hot 
planets (a < 0. 1 AU) are not shown in this figure. 



ticularly metal rich, even among exoplanet hosts, manifesting in 
an absence of Hot Jupiter host stars with a clear subsolar metal- 
licity. A number of Hot Jupiters around at least fairly low [Fe/H] 

0.2 to -0.4 hosts are now known, and statistical tests point to 

rather insignificant differences in the metallicity distribution of 
close-in planets and those on wider orbits (Ammler-von Eiff et 
al. 120091 1. On the other hand, planets of metal poor stars might 
have a tendency to have rather large semimajor axes outside 1.5 
AU (Santos et al.|2009). The basic problem is that our samples 
are still too small for definitive conclusions. 

The lack of a clear imprint of [Fe/H] on the semimajor axis 
might seem surprising at first glance from a theoretical point of 
view seen the significant influence of [Fe/H] on a start or Aa dis- 
cussed in the previous chapters. However, as we will demon- 
strate below, once these different effects are combined, one can 
actually understand such a weak dependance well on theoretical 
grounds. 



Figure 12 shows the distribution of final semimajor axes for 
planets more massive than 300 M m . The planets are binned into 
low, medium and high [Fe/H] cases. Note that each metallicity 
bin was normalized individually, so the absolute heights can- 
not be compared between different metallicity bins. The plot 
shows the nominal population (left panel), and a population with 
a; ce = 2.7 AU fixed (right panel). Note also that the "hot" plan- 
ets, i.e. those planets that reach the inner boarder of the compu- 
tational disk at about 0.1 AU are not included here. Due to the 
curvature of the feeding limit, also the bins out to a distance of 
about 0.2 AU are artificially lowered (see Fig.[Tji. The planets in 
the feeding limit are addressed below (Sect.[7]i. 

While the exact numerical values are different for the two 
populations, we see that the general shape of the semimajor axis 
distribution is similar, and consists first of a slow increase of the 
planet frequency with increasing semimajor axis for high [Fe/H] 
respectively an approximately flat part for lower [Fe/H], then an 



upturn to a maximum at a few AU (at 2-3 AU on the left, at ~ 1 
AU on the right), followed by a gradual decline at even larger 
radii. This is due to the fact that the core accretion timescale be- 
comes too long at large distances (Sect. 5.1.2 1. Our model does 



not include at the moment mechanisms like scattering that could 
bring planets to larger semimajor axes. This probably leads to an 
underestimate of the number of giant planets at large distances, 
which is important for direct imaging searches (e.g. Bonavita et 
al. 12009 l l. This issue will be addressed in forthcoming work. 



6.3.1. Imprint of the iceline 

For the nominal population, the location of the peak can roughly 
be understood by considering some relevant mean values (me- 
dian values give similar results) for the planets with a final mass 
> 300 M ffi : First, the mean Eo of these planets is about 480 g/cm 2 , 
corresponding to an a lcf , m 5.5 AU. As indicated by Fig. [8] the 
mean starting position should be somewhat outside of this dis- 
tance. The measured mean a stalt is indeed 6.5 AU. The mean 
disk lifetime of 3.8 Myr corresponds to a mean Aa w -3.8 AU 
(Fig. _ 



11 



so that we estimate a mean final semimajor axis of 2.7 
AU, which is the same as the mean value found in the simula- 
tions. This corresponds approximately to the peak of the blue 
line, which is the bin with the highest number of planets in it. 
Thus, the peak of the distribution is due to the fact that there is a 
typical locus from where planets come, combined with a typical 
distance over which they migrate. 

This imprint of the ice line into the semimajor axis distribu- 
tion is confirmed by the right panel, where ai Ce is at 2.7 AU for 
all disks. Compared to the nominal case, the upturn in the distri- 
bution is found at a smaller semimajor axis, namely at about 1 
AU which might be closer to the observed distribution and might 
also be sharper, at least for the high metallicity bin. These differ- 
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ences can be understood: The sharpness is a logical consequence 
because part of the spread in the typical a start > ai ce is now sim- 
ply eliminated by suppressing the dependence of ai ce on Sq. The 
fact that the upturn happens at a smaller distance simply comes 
from the fact that an a; ce = 2.7 AU is clearly less that the mean 
a; ce in the nominal case (5.5 AU). We conclude that there is a 
correlation between the position and shape of the upturn in the 
planetary semimajor axis distribution, and the thermodynamic 
properties in the protoplanetary disk, which is an interesting re- 
sult. 



Table 1. Mean values of the starting position a s tart> tire extent 
of migration Aa, and the final semimajor axis a for planets with 
a final mass larger than 300 M®, as a function of metallicity. 
Planets inside 0. 1 AU are excluded. 



Metallicity 


"start [AU] 


Aa [AU] 


a [AU] 


[Fe/H]> 0.2 


6.2 


-3.5 


2.7 


-0.2 <[Fe/H]< 0.2 


6.7 


-4.0 


2.7 


[Fe/H]< -0.2 


7.5 


-4.9 


2.6 



6.3.2. Comparison with observations and predicted decrease 

Note that we here only consider the condensation of ices as 
a mechanism causing a sweet spot for giant planet formation. 
Other (or additional) mechanism with comparable consequences 
could also be at work in the disk (cf. Schlaufman et al. 2009). 
One should also keep in mind that the temperature and solid 
surface density structure of protoplanetary disks is likely much 
more complex than assumed in the simple a models used here 
(Dzyurkevich et al. 12010b . This means that the result that the 
frequency of giant planets is an increasing function of the semi- 
major axis inside a few AU is a more robust prediction than the 
specific shape and cause of it. Observationally, after correction 
of the observational bias, an overall distribution of the semimajor 
axes of giant planets not dissimilar from the one here has been 
found recently (Mayor et al. 1201 11 1. No imprint of [Fe/H] is seen 
in the still small dataset, where no fine structures (like e.g. the 
small valley visible at low [Fe/H] just inside the upturn) could 
currently be seen. It should be noted that such fine structures are 
sensitive to specific model assumptions. 

The baseline of the RV surveys is also not yet long enough 
to determine if the theoretically predicted decrease at larger dis- 
tances (> 5 AU) exists. Such a decrease can in contrast be seen 
as a solid prediction of the core accretion theory. The observa- 
tional test of this will be an important constraint for the core 
accretion model, as it is linked to the core growth timescale. A 
complication arises from the fact that the gravitational interac- 
tion of several giant planets can lead to the scattering of bod- 
ies to larger semimajor axes. This evolutionary effect will blur 
the timescale limit set during formation. To still be able to use 
the observational constraint, formation and evolution through N- 
body interaction must be coupled self-consistently (Alibert et al. 
in prep.). 



6.3.3. Independence from [Fe/H] 

Comparing the three metallicity bins in the nominal case, we see 
that there are some visible, but not extreme differences. For the 
low metallicity bin, one sees that this distribution is rather fiat 
inside about 3 AU, then has a sudden upturn, followed again by 
a sharp downturn. Thus, the planets are mostly well confined to a 
particular region. For the high metallicity bin, and to a lesser ex- 
tent also for the intermediate metallicities, the increase at small 
semimajor axes is more gradual towards the peak, and also the 
decrease towards the even larger semimajor axis is less abrupt, 
which means that the high [Fe/H] distribution has more planets 
both at smaller and larger distances than the low metallicity one. 

These finding are due to the fact that giant planets around 
metal poor stars can only form within a well defined region, 
namely at the sweet spot for giant planet formation as was shown 
in Sect. 5.2 in contrast to higher [Fe/H], which can form all over 
the disk. Considering the location of the peak for the low metal- 



licity bin, one notes that the peak of giant planets around low 
metallicity stars is found somewhat further out by (about 1 to 3 
AU) than for the higher metallicity cases. This might be in ac- 
cordance with the observational findings of Santos et al. (|2009). 
This can also be understood in terms of our findings from the 
previous sections: At low metallicities, giant planets must come 
from outside the ice line (Sect. |5.2) , i.e. from a significant dis- 
tance. At the largest semimajor axes (a > 6 AU), the relative 
frequency of more metal rich planets is again higher. This is a 
consequence of the timescale effect: a lot of solids are necessary 
to form a core in time at large distances, and this is not the case 
at low [Fe/H]. 

These reasonings alone are however not yet sufficient to ex- 
plain the observed (and theoretical) semimajor axis distribution 
with its weak [Fe/H] dependence, as it would rather imply a 
strong absence of giant planets around low metallicity stars at 
small semimajor axes. The missing link is the anti-correlation 
of metallicity and the extent of migration (sect. 6. 1 1, due to the 
"braking effect" which stops migration quicker in metal rich 
disk, and the "collection effect" which enhances it in solid poor 
ones. Thus, while giant planets initially arise closer-in in high 
[Fe/H] environments, they migrate less, while for low [Fe/H], 
they start further out, but migrate more. 

These findings can be quantified (Table [TJ: For the nominal 
population, the mean a start is 6.2, 6.7 and 7.5 AU for the high, 
medium and low [Fe/H] group. But the mean Aa is -3.5, -4.0 
and -4.9 AU again for high, medium and low [Fe/H], which re- 
sults in a near cancelation of the two counteracting effects, so 
that the resulting mean final semimajor axis for the three metal- 
licity bins are all the same except for a tiny difference of 0.1 AU. 
Observationally, such three mean values appear indistinguish- 
able. Thus, a complex of interplay of migration and accretion 
explains the observed absence of a strong correlation of metal- 
licity and semimajor axis for extrasolar giant planets. 

7. Metallicity of close in Jovian and Neptunian 
planets 

7.1. Observed distribution 

Early detections made with high precision radial velocity ob- 
servations (Udry et al. (2006) suggested that the distribution of 
the host star metallicities of close-in lower mass Neptunian and 
Super-Earth planets follows a different pattern than observed 
for hot Jupiters, namely that their metallicity distribution is not 
shifted towards the higher metallicities. This result was sub- 
sequently discussed in several studies (e.g. Sousa et al. 2008; 
Ghezzi et al. I20 10l > and recently confirmed (Mayor et al. I201 It 
Sousa et al. l2011l . 

The latter authors find that the mean metallicity of FGK hosts 
stars is positively correlated with the mass of their most mas- 
sive planet. They find a mean [Fe/H] of -0.11, 0.04 and 0.10 dex 
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Fig. 13. Histogram and cumulative distri- 
bution of the metallicity of close-in syn- 
thetic planets. The upper two panels are for 
the nominal population (fi = 0.001), the 
lower two for a population with a faster 
0.6 type I migration (fi = 0.1). In all pan- 
els, green lines stand for massive planets 
(M > 100 M e ), the blue lines are planets 
with 100 > M/M e > 6 M e , and red lines 
are low mass planets (M/M e < 6), all in- 
side the feeding limit, i.e. inside about 0.1 
AU. The black lines serves for compari- 
son and shows the metallicity distribution 
of all synthetic planets (at all semimajor 
axes) and thus simply the distribution from 
which we draw the initial conditions. 
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for host stars having planets with masses between 0.01 and 0.1 
Ma , 0.1 and 1 M^ and > lM^, respectively. In a similar way 
show Ghezzi et al. (2010) that for FGK host stars, in the obser- 
vational data, there is an offset in the mean metallicity of +0.1 1 
dex by which Jovian planet hosts are more metal rich than stars 
which host Neptunian planets only. The probability that the two 
samples are drawn from the same parent distribution is however 
found to be 17%, i.e. clearly non-negligible. They give a mean 
metallicity for Neptunian planet only FGK hosts of 0.01 dex, and 
of 0.12 dex for Jovian planet FGK hosts. 



Table 2. Mean [Fe/H] of the initial conditions leading to differ- 
ent types of synthetic planets. 



Mean [Fe/H] 


Nominal 


/i = 0.1 


All planets (initial cond.) 


-0.02 


-0.02 


Hot, M/M e < 6 


0.00 


0.00 


Hot, 6 < M/M e < 100 


0.05 


0.06 


Hot, M/M s > 100 


0.08 


0.17 



7.2. Synthetic result 

It is obvious that this raises the question whether such a rela- 



tion exists for the synthetic population, too. From Sect. 4. 1 such 



a behavior is expected, but the problem is that there, the mass 
function is shown for all semimajor axes, while the observed 
planets with M < 30 M e are almost all inside of a few 0. 1 AU. 

the metallicity distributiorj^] of 



We therefore plot in Fig. 



13 



those synthetic planets that have migrated into the feeding limit 
at about 0.1 AU ("Hot" planets), for three different mass bins: 
M > 100Me (green lines), 100 > M/M e > 6 (blue lines) and 
M/Me < 6. The lower limit of 6 Mb approximately reflects the 
lower limit to which the model should yield masses which are 
not too much affected by the initial conditions (Paper I). Below 
this limit we can still conclude that there are additional low mass 



2 With metallicity of a planet we refer in this work always to the 
disk/stellar [Fe/H] of the initial conditions which lead to the formation 
of this planet, and not to the heavy element abundance in the planet 
itself. 



planets, but their final mass is not well described in the model. 
The upper two panels in Fig. 13 show the nominal population 
with a type I efficiency factor fi = 0.001, while the two lower 
panels show a non-nominal population with f\ — 0.1. Numerical 
values for the mean metallicity of the different sub-populations 
are listed in table [2] 

In all figures, the black line additionally shows for compari- 
son the [Fe/H] distribution of all synthetic planets, and thus the 
distribution from which the initial conditions are drawn. It is a 
Gaussian with a mean fi = -0.02 dex and a dispersion cr = 0.22 
dex (cf. Paper I). 

One notes that for both populations, the hot Jovian planets 
have a distribution that is more metal rich than the intermedi- 
ate mass and Neptunian planets or the Super-Earth planets. For 
the nominal distribution, the difference between the intermediate 
and the high mass bin is however extremely small, except for a 
possible excess of Hot Jupiters with a high metallicity > 0.2 dex. 
In terms of mean metallicities, the difference is just 0.03 dex. For 
the population with fi = 0.1, the situation is quantitatively quite 
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Fig. 14. Left panel: Detection probability as 
a function of metallicity and gas surface 
density. The plot shows the relative detec- 
tion probability as a function of [Fe/H] for 
disks which have a So meeting the condi- 
tion indicated in the plot, normalized by the 
f(jet([F e /H], all So) i- e - the detection proba- 
bility as a function of [Fe/H] alone, includ- 
ing all So. This quantity is plotted in the 
right panel. A curve in the left panel which 
increase towards the low [Fe/H] (blue lines) 
therefore means a less strong metallicity 
effect, while a curve that decreases to- 
wards the low [Fe/H] (red lines) indicates a 
stronger metallicity effect. The instrumen- 
tal accuracy is 10 m/s. 



different: Here, the mean metallicity of the Jovian mass bin is 
with 0.17 dex clearly higher (+0.11 dex) than the middle bin 
with the Neptunian planets (0.06 dex). For the lowest mass bin, 
which contains the proto-terrestrial planets, a vanishingly small 
offset to the distribution from which the initial conditions are 
drawn is seen, for both synthetic populations, as these planets 
have a mean [Fe/H] =0.00. 

We thus see, as expected, that the mean [Fe/H] increases 
with increasing planetary mass. Quantitative results are however 
dependent on poorly constrained model parameters like fi. The 
large difference for the mean metallicity of Jovian planets for the 
two populations stems in particular from the following: Giant 
planets forming in high [Fe/H] and low So environments inside 
the ice line do not usually migrate into the feeding limit pro- 
vided that type I migration is negligible (fi = 0.001), because of 
the very efficient braking effect in type II migration at such small 
semimajor axes (sect. 



6.1 



. They rather stay at distances between 
roughly 0.4 and 1 AU. For fi =0.1 this is different: there, type 
I migration, which has a migration rate that increases in con- 
trast linearly with planet mass, brings the quickly growing cores 
starting inside the ice line already so close to 0. 1 AU, that many 
planets forming in such a high [Fe/H] environment eventually 
end up as Hot Jupiters, causing thus the high mean metallicity. 

For a third population with fi = 0.001, but only partially 
suppressed type II migration, one finds results for the metallici- 
ties that approximately lie between the two cases discussed here. 
This means that the metallicity distribution of Hot Jupiters is a 
measure of the efficiency of both type I and type II migration. In 
summary we see that the synthetic populations reproduces the 
general observed trend, but that the more specific results depend 
on uncertain model settings like the efficiency of migration. 



8. Planet frequency 

Comparing observed and theoretical detection probabilities Pdet 
of planets as a function of some input (disk) variable is a clas- 
sical application of population synthesis calculations (e.g. Ida 
& Lin l2004bl Korn et et al. 120051 Dodson-Robinson et al. l2006l 
Matsuo et al. l2007t . Here we study the frequency of giant planets 
as a function of [Fe/H], Mdi s k and Tdi s k- 



8.1. [Fe/H]: Influence of the gas mass 

In Paper II, we compared the synthetic "metallicity effect" (the 
increase of the detection rate of giant planets with host star 



[Fe/H]) to the observed one and found that the synthesis can re- 
produce this important observational constraint. Here we focus 
on another aspect, namely how the "metallicity effect" depends 
on the amount of gas present in a protoplanetary disk. 

For simplicity, we here again assume that a planet can be 
detected by the radial velocity method if it induces a velocity 
semi-amplitude K larger than 10 m/s, and if its orbital period is 
less than 10 years. We thus focus on giant planets. 

left panel, shows the detection probability as a 



Figure 
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function of [Fe/H] if additionally the gas mass S of the corre- 
sponding disk is either larger than some lower limit (blue dashed 
lines), or smaller than some upper limit (red dotted lines). In all 
cases, the detection probability Pd et ([Fe/H],So) obtained in this 
way was divided by the detection probability at a given [Fe/H] 
but including all So, i.e. the (usual) detection probability as a 
function of [Fe/H] alone. This quantity is shown in the right 
panel of Fig. [14] Like that, the differential influence of the gas 
mass on P^et is shown. 

One notes that for disks which have a So above a certain 
threshold (blue lines), i.e. where the most gas poor disks are ex- 
cluded, the normalized probability is larger than 1, which means 
that more planets are detected. When only a moderate lower limit 
is used (So > 100 g/cm 2 ), the increase is only small and nearly 
identical for all [Fe/H]. However, the higher the So threshold be- 
comes, the more a dependence on [Fe/H] develops: The relative 
increase of Pj et is highest for the subsolar metallicities. In other 
words, the metallicity effect is weakened, if the protoplanetary 
disks are gas rich. A mirrored effect is seen if only disks with 
a mass smaller than some S (red lines) are considered, so that 
the most gas rich disks are now excluded: Here, the relative de- 
crease of fdet is strongest for the low [Fe/H] bins, which means 
that the metallicity is strengthened if the protoplanetary disks are 
gas poor. 

These findings are clearly a consequence of the "compen- 
sation" effect discussed theoretically in Sect. 5.3 In a gas rich 



environment, a lower [Fe/H] is sufficient to allow giant planet 
formation. 



8.2. Disk gas mass 



Figure 15 shows the percentage of stars (initial conditions) 



where the formation of a giant planet with a mass of at least 
100M ffi is possible, as a function of the initial disk gas mass 
measured as [Md/Msn]- The lowest value of about -0.6 (So = 50 
g/cm 2 ) corresponds to an initial content of gas within the com- 
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Fig. 15. Percentage of stars having a giant planet with a mass of 
at least 100 M e as a function of the initial disk gas mass. The 
dotted line scales as M^ k . 



putational disk which of about 0.004 M Q . The upper limit of 
about 0.7 (Eq = 1000 g/cm 2 ) corresponds to a disk with an ini- 
tial mass of 0.09 M e . As the individual primordial disk mass 
is unknown for the Gyr old stars around which extrasolar plan- 
ets are typically discovered, we directly plot the real frequency, 
without taking into account any observational bias, in contrast to 
the results for [Fe/H] . 



As expected from sections 4.1 one sees that the frequency 
of giant planet is increasing with [M D /Ms N ]. For [M d /Msn] ^ 
-0.1, the increase is found to be slightly stronger than linear 
with the disk mass, as indicated by the dotted fitting line which 
correspond to a fraction of stars with giant planets given as 
O.O4x(Mdi s k/O.O17M ) L2 . For lower disk masses, this fit overes- 
timates the number of giant planets. We note that two quantities 
which are important for giant planet formation scale with sim- 
ilar exponents: the core accretion rate scales linearly with the 
planetesimal surface density Ss which is itself proportional to 
the disk gas mass (Ss = /d/g^oX while the isolation mass is 
proportional to Eg 5 (e.g. Lissauer et al. H993t . 

For the most massive disks considered in this work, one finds 
a frequency of about 28%. This is in fair agreement with the re- 
sult of ~ 35% of Dodson-Robinson et al. (2006) at the same 
disk mass. At lower disk masses we find somewhat higher fre- 
quencies than these authors. The reason for this difference could 
be that in Dodson-Robinson et al. (2006), migration was not in- 
cluded, so that higher surface densities are necessary (cf. sect. 

HD- 

8.3. Disk lifetime 



Figure 16 shows the percentage of stars with a giant planet 
(M > 100M ffl ) as a function of disk lifetime Tdi s k- The solid 
line shows the dependence for the nominal population. Not sur- 
prisingly, there is a positive correlation of the disk lifetime and 
giant planet occurrence, as in longer living disk, cores have more 
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Fig. 16. Percentage of stars having a giant planet with a mass 
of at least 100 M e as a function of the disk lifetime Tdisk- The 
solid line shows the nominal population, while the dotted line 
corresponds to a population where disk mass and disk lifetime 
are not correlated (cf. sect. |9jl. The dashed lines scales as 



r disk' 



time to grow to a size where they can trigger gas runaway accre- 
tion. The dashed lines shows an approximate fit to the frequency, 
given as 0.075 x (Tdisk/3 Myr) 2 . The strong dependence on Tdisk 
is caused by the double role of Tdisk both as a threshold, and as a 
factor directly influencing the planet mass, as described in Sect. 

For the disk with the longest Tdisk ~ 6.5 Myrs, the percentage 
of stars with a giant planet is remarkably high, namely about 
40%. For comparison, Dodson-Robinson et al. (2006 ) find about 
30% for such a Tdisk- This high fraction is however dependent 
on a disk mass-disk lifetime correlation inherent in the initial 
conditions, as addressed in Section|9]below. 

Recently, it was discussed whether extrasolar planet host 
stars show a stronger depletion in lithium compared to non- 
planet hosts. No general agreement exists currently on this 
question on observational bases, with analyses both in favor 
(Israelian et al., 120091 Sousa et al. 120101) or against (Baumann 
et al. 120101 Ghezzi et al. 120101) such a correlation. On theo- 
retical bases, several mechanism can enhance lithium depletion 
(e.g. Baraffe & Chabrier 2010). A possible explanation based 
on the star rotational history which is interesting in the context 
here has been put forward by Bouvier (2008): In this model it 



is assumed that long disk lifetimes lead via disc locking to a 
slower initial rotation of the star at the ZAMS. This causes a 
phase of stronger differential rotation between the radiative core 
and the convective envelope, which is in turn responsible for an 
enhanced lithium depletion. 

If the correlation of enhanced lithium depletion and the pres- 
ence of massive planets can indeed be confirmed, then the lines 
of reasoning of Bouvier (2008) provide a way of relating disk 
lifetime and the formation probability of giant planets. The result 
shown here that giant planets preferentially form around disks 
with a long lifetime fits well in this picture. Note that in the di- 
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rect gravitational collapse model for giant planet formation, this 
is not the case. 

The abundance of lithium would then provide a link which 
is observable today between primordial disk properties (Tdisk) 
and planet occurrence in similar, although much more indirect 
way than [Fe/H]. We comment that long disk lifetimes are in 
particular needed for giant planet formation when the metallicity 
is on the low side ([Fe/H]< -0.2), while at high [Fe/H], short 
lifetimes (Tdisk 55 2 Myrs) can be sufficient, as demonstrated in 
Section|5 .3.2| so that the planet frequency-Td; s k correlation also 



depends on [Fe/H] . 



9. Disk mass - disk lifetime correlation 

When generating initial conditions, we first draw in the nominal 
procedure the gas surface density S and then, independently of 
it, the photoevaporation rate M w . This procedure has the side 
effect that disks with a high initial gas surface density So will 
on average also have a higher lifetime r^k, as for the same M w 
more massive disks life longer. Thus, there is a positive correla- 
tion between disk mass and lifetime. 

To assess the consequences of this correlation, we have gen- 
erated an alternative set of non-nominal initial conditions in the 
following way: After drawing So, we directly draw a r^k from 
a specified distribution and then choose the M w that is necessary 
to get this Tdisk, for the given So (and the fixed a). Like that, 
no correlations exist any more between So and r^k- Instead, So 
and M w are now positively correlated, which could arise in a sit- 
uation where photoevaporation by the central star is dominant, 
even though that then, strictly speaking, the functional form of 
the photoevaporation would be different than assumed in our 
disk model (Matsuhama et al. !2003b . 

The distributions we have used for Tdisk for the non-nominal 
case are either a uniform distribution of lifetimes between and 
6 Myrs which corresponds to the linearly decreasing fraction of 
stars which have a JHKL excess as a function of mean cluster 
age in Haisch et al. (2001 ), or an exponential decrease of this 
fraction with time as exp(-f/2.5 Myr), which is based on more 
recent compilations of disk lifetimes as in Mamajek (2009) or 
Fedele et al. d20T0l . 

9. 1 . Effect on the giant planet frequency 

We have found that the impact of using this different prescription 
only has minor impact on the general properties of the synthetic 
population as compared to the nominal case. This regards in par- 
ticular the distribution in the mass-distance plane, the migration 
behavior or the general shape of the PIMF. The difference be- 
tween the linear and the exponential distribution for Tdisk is even 
smaller. 

Only when we directly plot the frequency of planets as a 
function of disk lifetime as in Fig. 
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a significant difference 
becomes visible. In this figure, the percentage of stars with a 
giant planet for the non-nominal case with the linear distribu- 
tion is shown as the dotted line. One notes that while for the 
nominal case (solid line), a clear increase of the giant planet fre- 
quency with disk lifetime is seen for Tdisk ^ 4 Myrs, this is not 
the case for the population without a correlation between Tdisk 
and Mdisk- The reason for this is that in the non-nominal case, 
there are many disks with a long lifetime, but still only such a 
small So that giant planets cannot form, which reduces the frac- 
tion. Before, disks with a high Td; s k usually also had a high So, 
which contribute both to successful giant planet formation. 



10. Summary and conclusion 

We have extensively studied for solar-like stars the influence of 
some of the most important properties of their protoplanetary 
disks (disk metallicity [Fe/H], disk (gas) mass Mdisk and lifetime 
Tdisk) on the planets growing in them. We have found numerous 
correlations which we summarize here. 

We first give a list of the most central correlations between 
disk and planetary properties as found within the core accre- 
tion paradigm, which depend less sensitively on model settings. 
From an observational point of view, we have to distinguish two 
classes: First correlations which can be tested against observa- 
tions, as they involve the (stellar) [Fe/H]. Second correlations 
which are linked to the disk mass and lifetime. These primor- 
dial disk properties cannot be directly observed, even though that 
some pathways to at least approximately determine them might 
arise in future. The correlations involving the metallicity are: 

1 . The planetary initial mass function for metal rich disks con- 
tains a higher number of giant planets, but not of a signifi- 
cantly higher mass. Metallicity mostly acts as a threshold for 
giant planet formation, but is not correlated with the mass of 
giant planets, except for rare special cases (pointQ. 

2. No clear metallicity effect (positive correlation of the de- 
tection rate of planets with stellar [Fe/H]) is found in the 
Neptunian mass domain. At even lower masses, an anti- 
correlaton could exist. 

3. Observationally, at a low radial velocity precision i.e. for 
giant planets, no clear correlation of [Fe/H] and the plane- 
tary mass distribution is visible. At high precision (~ 1 m/s) 
in contrast, a systematic anti-correlation of [Fe/H] and the 
Neptunian to Jovian planet ratio becomes detectable: The 
lower [Fe/H], the higher the ratio of Neptunian to Jovian 
planets. 

4. The most massive companions (> 10 - 20 M^) cannot form 
by core accretion at clearly subsolar [Fe/H], as core growth 
takes then so long that disks are significantly depleted once 
the cores trigger gas runaway accretion. In the supersolar 
regime, maximal planet masses are in contrast independent 
of [Fe/H]. Such planets are however anyway very rare. 

5. The absence of a strong correlation of metallicity and the 
semimajor axis distribution of giant planet is explained like 
this: Around metal poor stars, giant planet cores can only 
form in gas rich disk, and at large distances a s tart <: ^ice- 
These high gas disk masses cause migration over larger dis- 
tances, tending to cancel out the large a s tart- In metal rich 
disk it is the contrary: giant planets can form also at smaller 
distances. But as low gas masses are sufficient, planets mi- 
grate also less. These correlations of accretion and migra- 
tion explain the observed weak dependence. In short: at low 
[Fe/H], planets start further out, but migrate more, while at 
high [Fe/H] they start closer in, but migrate less. 

The correlations involving the disk mass and lifetime are: 

1. The planetary initial mass function for massive gas disks 
contains more giants planets of a higher mass, but less giants 
of a lower mass. Disk gas masses and giant planet masses are 
correlated. 

2. Disk lifetimes act both as a threshold (like [Fe/H]) for giant 
planet formation but also influence the final total mass (like 
Mdisk), so that they have a strong influence on the giant planet 
population. 

3. Maximal planetary masses and disk mass are in contrast to 
the situation for [Fe/H] correlated over the full range of disk 
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masses. Except for very low disk masses, maximal possible 
planet mass and disk mass correlate linearly. There is a wide 
spread in the efficiencies of converting disk gas in envelope 
mass, but to order of magnitude it is 10%. 

4. Low metallicities (i.e. dust-to-gas ratios) can be compen- 
sated by high disk masses to allow giant planet formation, 
and vice versa ("compensation effect"). The lowest neces- 
sary total solid content in a disk for giant planet formation is 
about 0.4 to 0.5 M^. 

5. The metallicity effect is weaker when disk gas masses are 
high, and stronger if they are small. This is a consequence of 
the metallicity-disk gas mass "compensation effect". 

6. Near the solid surface density threshold necessary for giant 
planet formation, long disk lifetimes of at least 3-4 Myrs are 
necessary. This could indicate that the Solar Nebula probably 
had a rather long lifetime. 

We have also found the following correlations, which are ei- 
ther of second order, or could be dependent in a more significant 
way on specific model assumptions: 

1. There is an imprint of the disk mass distribution on the upper 
end of the planetary mass function. 

2. The chain of correlations of low [Fe/H]=> high necessary gas 
mass => large a; ce => large typical a stalt means that in low 
[Fe/H] disks, embryos of giant planets come typically from 
further out. The opposite applies for high [Fe/H]. 

3. For giant planet formation, low metallicities can not be com- 
pensated by high gas masses ad infinitum, at least if higher 
mass disks have an ice line further out due to stronger vis- 
cous dissipation: The sweet spot for giant planet formation 
at a s tart~<2ice moves then so far out that the increasing core 
growth timescale can no more be compensated for by in- 
creasing the disk mass. 

4. The semimajor axis distribution of giant planets outside 0. 1 
AU first consists of a nearly fiat part, followed by an upturn. 
The location of the upturn probably depends on the location 
of the ice line. It is due to the fact that there is a typical 
locus from where giant planets come, combined with a typi- 
cal distance over which they migrate. Disk thermodynamical 
properties thus leave an imprint in the planetary semimajor 
axis distribution. 

5. The extent of migration of giant planets Aa = a-a s tart is large 

( 6 AU) in gas rich disks with a low [Fe/H]: Torques are 

strong, a start are large, and cores must migrate over a large 
distance until they have collected sufficient planetesimals to 
trigger gas runaway accretion and slow down in slow planet 
dominated type II migration. The extent is minimal (< -1 
AU) in gas and solid rich disk and if the embryo starts in- 
side fli ce : giant planets grow rapidly and local disk masses are 
small, so that planets soon stop due to their inertia. For low 
mass disks Aa — 3 AU, and for disks similar to the MMSN, 
Aa —A AU. 

6. The extent of migration of giant planets increases with disk 
lifetime Tdisk ■ Very roughly, there is a linear correlation. For 
disks with the shortest * 0.5 Myrs still allowing giant 
planet formation, Aa is — 2 AU, while in the longest living 
disks (T^sk ~ 7 Myrs), Aa — 5.5 AU. There are important 
systematic departures from this correlation depending on the 
solid content of the disk. 

7. Observationally, Hot Jupiters are found around higher metal- 
licity hosts than lower mass close-in Neptunian and Super 
Earth planets. This is also found in the synthetic population. 
The degree of difference correlates positively with the effi- 



ciency of type I migration and varies between 0.03 and 0. 1 1 
dex for the mean [Fe/H] values. 

Finally, with the following points we mostly confirm ear- 
lier results concerning the relationship of disk properties and 
planet properties. The quantitative results often differ to a cer- 
tain degree, but qualitatively they are the same (cf. Ida & Lin 
l2004all2004Tl Dod son-Robinson et al.l2006l Kornet et al.[2006; 
Matsuo et al. l2007l Thommes et al. 120081 

1. Comparison with observation shows that the very large ma- 
jority of known extrasolar planets lies in the [Fe/H] -mass do- 
main populated by the core accretion formation mechanism. 

2. The sweet spot for planet formation is at about 5 AU. The 
minimal necessary surface density of planetesimals to allow 
giant planet formation is about Ss = 6 g/cm 2 there. Both in- 
and outside this distance, higher are necessary. Inside, to 
overcome the small isolation mass, outside to overcome the 
long core growth timescale. 

3. At metallicities [Fe/H]< —0.1, giant planets can only form in 
an annulus of 2-3 AU in width outside the ice line a; ce . With 
increasing [Fe/H], giant planets can also form in- and clearly 
outside a; ce , but the ice line remains the preferred starting 
position. Only at the highest [Fe/H]> 0.4, the location of aj ce 
becomes less important and giant planets can form all over 
the disk from ~ 1 to ~ 20 AU. 

4. The fraction of stars with a giant planet as a function 
of initial disk mass is found to be approximately 0.04 x 
(M dlsk /0.017M o ) 12 . 

5. The fraction of stars with a giant planet as a function of disk 
lifetime is found to be roughly 0.075 x (Tdi s k/3 Myr) 2 . A pos- 
sible link between giant planet occurrence and lithium deple- 
tion in the host star could be an observational manifestation 
of this correlation. The influence of Tdisk on the frequency of 
giant planets is however much reduced if the disk mass and 
lifetime are not positively correlated, as implicitly assumed 
in the nominal population. 

This long list of correlations clearly indicates the strong de- 
pendence of planet formation on protoplanetary disk properties. 
It calls for profound improvements of the disk models used in 
planet formation simulations, for a rigid inclusion of observa- 
tional results on disks in the initial conditions, and for the incor- 
poration of observations of disks like SEDs as additional con- 
straints for planet formation models. 

Such improvements of the initial and boundary conditions 
for planet formation will in the end allow to better understand 
the planet formation process itself. 

Note: The numerical data of the nominal syn- 
thetic population can be obtained in electronic form at 
http://www.mpia.de/homes/mordasini/Site7.html 
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